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THEORETICAL AND EXPERIMENTAL STUDIES OF REENTRY PLASMAS 


By Michael G. Dunn and Sang -Wook Kang 
Cornell Aeronautical Laboratory, Inc. 

SUMMARY 

A viscous shock-layer analysis has been developed, and applied to the 
calculation of nonequilibrium -flow species distributions in the plasma layer 
of the RAM vehicle at high altitudes. The theoretical electron-density results 
obtained are in good agreement with those measured in flight. The flight 
measurements were obtained using electrostatic probes that protruded well 
into the shock layer. In addition, the theoretically obtained heavy-particle 
translational temperatures appear to agree fairly well with the electron tem- 
peratures that were measured in the flight experiments using voltage- swept thin- 
wire electrostatic probes. The influence of the reaction rate coefficients on 
the calculated electron densities has been assessed and shown to be within the 
uncertainty in the flight data. The theoretical results demonstrate the importance 
of including in the chemical model the positive ions N£ , » N , and O , 

in addition to NO^, for the high altitudes and velocities considered here. 

A circular- ape rture flush-mounted antenna, for which a theoretical 
description of the antenna-plasma interaction has been developed, was used 
to obtain a comparison between theoretical and experimental antenna admittance 
in the presence of ionized boundary layers of low collision frequency. The 
antenna was located in a flat-plate model which was in turn located in the 
expanding-flow environment of a r eflected- shock tunnel. The electron-tem- 
perature and electron-density distributions in the plate boundary layer at the 
antenna location were independently measured using voltage- swept thin- wire 
Langmuir probes for one of the test conditions. The antenna admittance 
was measured using a four-probe microwave reflectometer and these 
measured values were found to be in good agreement with those predicted 
from the theory. The response of this antenna was found to be very 
sensitive to changes in the profile and magnitude of electron density for 
0 . 1 <n e / n e ) cr < 1.2. Beyond n e /n e ) cr of 1 . 2 the phase rotation continued 
to provide a measurement of electron density. Measurements were also 
performed with another type of circular-aperture antenna which did not 
satisfy all of the constraints of the theoretical model, yet good agreement 
was obtained between the calculations and the experimental results. 


Voltage-swept thin-wire electrostatic probes have been used to obtain 
the electron density in the free- stream flow of a ref lected- shock tunnel, and in 
the boundary layer over a flat plate located in this expanding-flow environment. 
The voltage sweep was sufficient to permit electron-density determination from 
both the ion- and electron-current portions of the probe characteristic while 
the electron temperature was determined from the electron- retarding region. 

The electron density determined from these two portions of the probe character- 
istic were in good agreement with each other and with independent microwave 
measurements for the free-stream experiments. Laframboise's free-molecular 
flow theory was used to deduce the number-density results presented here. 

A theoretical analysis has been completed which permits the calculation 
~>f the nonequilibrium, viscous shock-layer flow field for a sphere-cone body 
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at high altitudes. Results are presented for two different bodies at several 
different altitudes illustrating the influences of bluntness and chemical non- 
equilibrium on several gasdynamic parameters of interest. Plane-wave trans- 
mission coefficients were calculated for an approximate space- shuttle body- 
using a typical trajectory. These approximate results suggest that for fre- 
quencies less than 2200 MHz, transmission difficulties will probably be encoun- 
tered between 200 and 300 kft on the basis of pure air ionization only. The 
additional influence of ablation product ionization was not assessed. 

Some of the thin-wire electrostatic probe data obtained after the probes 
were retracted into the base region were analyzed and electron temperatures 
and electron densities were deduced over an altitude range of approximately 
190, 000 to 180, 000 feet. The results appear to be reasonable on the basis of 
electron densities and heavy-particle translational temperatures computed at 
the vehicle shoulder just prior to base-region expansion. 
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1. INTRODUCTION 


This document is the final report for Contract NAS 1-10674 which 
was initiated in April 1971. The work requirement of this contract consisted 
of five distinct problem areas each of which is sufficiently self contained 
that the results obtained can be reported separately. 

In the remainder of this report, the comparison between the theoret- 
ical and measured electron-density distributions for the RAM vehicle at 
high altitudes will be presented in Section 2. In Section 3, the experimental 
and theoretical studies performed using the RAM S-band atennas to determine 
electron density are described. The results of an experimental program 
using positively biased electrostatic probes to obtain electron density m 
collisionless flows are discussed in Section 4. A theoretical analysis which 
permits the calculation of the nonequilibrium, viscous shock-layer flow field 
for a sphere- cone body at high altitudes is presented in Section 5. This 
discussion includes calculated plane-wave transmission coefficients for an 
approximate space- shuttle body for frequencies from 10, 000 to 100 m z. n 
Section 6 the results of the post- retraction electrostatic-probe results are 
presented. The report concludes with Section 7 in which our suggestions for 
obtaining a description of the high altitude, three-dimensional space- shuttle 
flow field using the relatively simple integral method are presented. It was 
not within the scope of this contract to complete this particular study. The 
purpose was only to formulate the problem which was completed. 


2 THEORETICAL AND MEASURED ELECTRON-DENSITY DISTRIBUTIONS 

FOR THE RAM VEHICLE AT HIGH ALTITUDES 


2. 1 Introduction 

For hypersonic flight vehicles at high altitudes, the flow field sur- 
rounding the body, i. e. , the shock layer, is fully viscous, with the result 
that the conventional separation of the inviscid layer and the thin boundary 
layer is not warranted. For example, for a vehicle with a 12-inch nose 
diameter, the thin boundary-layer assumption breaks down at altitudes above 
200, 000 feet and a separate analysis is needed which accounts for the thick- 
ened shock-transition zone and the fully viscous fluid within the shock layer 
along with the effects of chemical nonequilibrium. Analyses are available 
for the nose region of a blunt vehicle (Refs. 1-4). It is the purpose of this paper 
1) to analyze theoretically the chemical nonequilibrium, ionized flow field 
surrounding a blunt-nosed entry vehicle at high altitudes, and 2) to compare 
the theoretical results with available in-flight measurements of electron- 
density distributions in the plasma layer. 
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The present analysis represents an extension of the basic treatment 
developed in Ref. 2 to the case of a body composed of a spherical nose region 
and a conical afterbody section for which the radius of curvature is infinite. 

For quantitative results, the extended theory has been applied to a NASA- 
developed RAM-C vehicle (hemisphere - 9° cone) for which large amounts 
of actual flight data are readily available (Refs. 5-12). For low altitudes, these 
measurements compared well with the theory developed (Refs. 13, 14) from a 
boundary-layer point of view. At high altitudes, the viscous-layer theory developed 
in the present paper will be shown to compare favorably with the flight data 
in terms of the electron-density profiles across the viscous shock layer at 
the base of the sphere-cone vehicle. 

The flight data reported here were obtained as part of the Langley 
Research Center's Project RAM (Radio Attenuation Measurements) which has 
for several years conducted in-flight measurements in order to improve under- 
standing of the causes and methods of predicting (Refs. 13-14) or preventing 
(Refs. 5-6) the radio- signal degradation during earth entry of a space vehicle. 
Electrostatis probes projecting from the spacecraft into the plasma layer have 
been flown in an attempt to measure in-flight electron-density distributions 
(Refs. 7,11). On-board diagnostic antennas (Ref. 8), reflectometers (Ref. 9). 
and VHF antennas (Ref. 10) have also been used to obtain information about the 
surrounding plasma. 

The theoretical analysis will be discussed in Section 2 and followed in 
Section 3 by a discussion of the chemical model and associated reaction rate 
coefficients. The experiment and diagnostic techniques are discussed briefly 
in Section 4. In Section 5, the theoretical results are compared with the 
flight data and the influences of uncertainties in the cheinical model and reac- 
tion rate coefficients are assessed. 


2. 2 Theoretical Analysis 

In order to present the analysis in a self-contained and coherent 
manner, the major aspects involved in the development of the theoretical 
analysis will be included here. More complete details may be found in 
Ref. 2. 

The assumptions made in the theoretical analysis are 1) a thin shock 
layer, 2) two-layer model of Cheng (Ref. 15), 3) constant Prandtl and Schmidt 
numbers, 4) binary diffusion due to concentration gradient only, 5) negligible 
changes in the flow properties due to the chemical reactions in the flow field, 

6) ambipolar diffusion for the electrons and ions. These assumptions have 
been made wholly or in part in the previous analyses and the justification for 
these assumptions can be found in Refs. 15-19- 

A. Basic Equations 

Based on the above assumptions the governing equations for the 
viscous, ionized shock layer in the merged-layer regime become (Refs. 2, 15) 
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These equations express the conservation of mass, momentum, energy, and 
chemical species, respectively, along with an equation of state. (See Fig. 1 
for a description of the flow field. ) The eleven chemical species considered 
are: O , N ' O , N , NO , NO + , O , N\ Ch , N ? ’and e . With 

the two-layer model of Cheng (Ref. 15), these equations are applied to the viscous 
shock layer allowing for the diffusion of the species into the shock-transition 
zone by using the modified Rankine-Hugoniot conditions. Analysis of the 
entire flow field performed in Ref. 17 reinforces the applicability of the two- 
layer model, and only the viscous shock layer will be considered here. 

The boundary conditions for a fully catalytic solid wall are (Refs. 2, 15): 
at = o j a - 0 = a j H — H b (X) , pp ~ jo^CX) 

( 0. , , N z ) (7) 
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The molecular concentrations are expressed in terms of the other species 
by locally conserving the elemental composition. In the case of pure air of 
present interest, this yields 
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where C N ( oo ) is taken to be 0. 767 and C 0 Coo) = 0. 233. 
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B. Chemical-Reaction Model 


The chemical reactions used simultaneously for dissociation and 
ionization are given and discussed in detail along with a table in Section 3. 
It may be mentioned here that eleven reactions are used involving only 
neutral species and fifteen reactions for charged species. 

C. Method of Analysis 


In analyzing the Eqs. (1-6) with the boundary conditions and the 
chemical reactions, the physical coordinates are transformed as follows : 
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In addition^ stream function \jr is introduced such that dy^/d% = - (1 + j)(TC r) Qlf, 
and 9^"/ dij. — (1 + (TC r y Q uu , which satisfies the continuity equation (1). 

A dimensionless stream function f ( £ , F ) is now obtained by putting 

yr = (1 f j ) o CTZ r )"* f , which yields the relationship fl-P/SF = U& CW/S/^ 

where U = U- / ( coo /S ) . 


Transformation of Eqs. (2-5) yields, after some rearrangement: 
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where Sc,; = ? t h •» H b /^ ao , & = ( /V - ^^/(^oo •> 

and K 2 is the rarefaction parameter due to Cheng (Ref. 15). The term ^ is a 
measure of the thickness of the viscous shock layer in the transformed plane. 

The Eqs . (11) - (14) should be considered simultaneously with 
appropriate boundary conditions. However, the previous results obtained, 
especially those by Chung, Holt and Liu, for the stagnation merged layer, 
indicate that the flow-field properties such as the pressure, normal velocity 
and enthalpy change very little despite the presence of chemical reactions m 
the flow. In the present study, the fluid dynamics equations, Eqs. (11) - y 5 ’’ 
are decoupled from the species conservation equations, Eq. (14), in an effort 
to simplify the analysis. Thus the flow properties such as the streamwise 
velocity, the pressure and the total enthalpy will be first obtained from 
Eqs. (11) - (13), so that they may be used as inputs to the species conservation 
equations. The static temperature, on the other hand, is sensitive to the 
chemical reactions present in the viscous shock layer, and must be considered 
along with the species conservation equations. 
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Hence, the problem now is to solve simultaneously the species- 
conservation Eq. (14) for O , N , NO , NO + , 0 + , N + , O 2 and N 2 + 
The electron-number density level may be obtained from the relationship 


n 


e * ^ n ^ {jt = ionized species), which gives in particles per c. c. : 


n = <? 


m, 


(15) 


where N e is Avogadro's number. 


In order to obtain solutions to the present problem, the Karman- 
Pohlhausen integral-method approach has been used because of the relative ease 
of application without unduly sacrificing the accuracy of the results (Refs. 2, 15, 
17-20). In addition, the method has been shown to be applicable to the case of 
nonequilibrium flow (Ref. 21), and to the case of merged-layer flow (Refs. 2,20). 
Thus, integrating the species conservation Eqs. (14) from F = 0 to F = 1 , we 
obtain after a series of rearrangements, 
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where * 2 cm / Of UC^dF and £. = (dC^/dF)^, 

and i = O, N, NO, °NO+, O , N , and . 

The various species profiles and the static temperature profile 
to be used in the above equations are, for a solid wall (see Ref. 2 for details 
of derivation): 
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By substitution of the Eqs. (17) and (18) in the definition of -fl A 
along with the other flow properties obtained from Eqs. (11) - (13), we 
obtain expressions for these quantities in terms of the unknown parameters 
£-. Other terms such as K , t b , G , etc. , ai;e specified as 

known input quantities. Since the species production term contains 

not only the temperature term but other species terms as well, it is neces- 
sary to solve Eqs. ( 1 6) simultaneously. Thus the problem is to determine 
from Eqs. (16) the unknown parameters E A as functions of the streamw 
distance £ . The term E 2 denotes the local concentration- gradient 

parameter in the transformed coordinate system. Solutions were obtained 
by integrating Eqs. (16) on a digital computer along the body from the stagna- 
tion region through the sphere-cone junction to the afterbody conical section 
by using the Adams-Molton predictor-corrector method. In order to simpli- 
fy the analysis while accounting for the change in shock shapes along the 
body, the distribution of the shock angle was initially assumed and updated 
from the results. Based on the cases considered in the present analysis, two 
iterations were usually sufficient to obtain correct shock shape along the body. 
The solutions are obtained in terms of the flow variables such as the veloci y 
and temperature, and in terms of the chemical species distributions m the 
shock layer and along the body surface. 


In order to test the validity of the theory, the present analysis was 
applied to a RAM-C vehicle (9° sphere-cone) because of the availability of 
the actual flight data on the vehicle, especially in terms of the electron- 
concentration distributions at high altitudes where the flow is m the viscous 
shock-layer regime. Solutions are obtained at several altitudes for compari- 
son between theory and experiment. This will be discus sed more fully m 
Section 5. 


2. 3 Chemical Model Used in Analysis 

The initial chemical model used in this study was composed of 
12 species and 64 chemical reactions. Subsequent calculations performed 
for the altitude and velocity regime of interest here made it possible to 
eliminate many of the unimportant reactions from the model. The reaction 
rate coefficients were perturbed by factors felt to be consist n ___ Hnns \ 

uncertainty (10 to 100 for many of the ion-atom and ion-molecule reactio » ) 

Drior to elimination of a particular reaction. The final chemical model u 
in the analysis is given in Table 1. The model included the neutral species 
N 2 , N, O z , O, and NO and the charged species N 2 +, N , 0 2 , O , NO • 
and e- . For the high altitude conditions of interest here (above 230 , 000 ft), 
calculations performed at NASA indicated that ablation-product ionization 
should not be a problem . 

The reaction rate coefficients for the reactions involving only neutral 
species (reactions 1-6 and 22-26) were taken from ^ + ^ & 

coefficients for the deionization of NO+, 0+, N , O 2 and N 2 
(reactions 7-10 and 13) were taken from the experimental results reported 

in Ref s . 29 - 32 . 
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Reaction rate coefficients for the charge-transfer and ion-molecule 
reactions (reactions 11* 12, 16-20) are not nearly as well known as those for 
the reactions discussed above. This lack of information is unfortunate be- 
cause these reactions do play an important part in the nonequilibrium process 
and can have a significant influence on the relative species distributions. 
However, for the calculations reported here, the only reactions of this class 
that were important. were reactions 11, 12, and 19* Of these, the rate co- 
efficients for reactions 11 and 12 were used in the data correlation of Refs. 31 
and 32 and are felt to be reasonably well known. Perturbation of the rate co- 
efficient of reaction 19 by a factor of 10 will change the relative concentrations 
of 02"*" and NO* somewhat but it will not influence the electron density signifi- 
cantly because of the fast two-body dissociative recombination of both NO* and 
0 2 * • Therefore, even though there is a relatively large uncertainty associated 
with the reaction rate coefficient for the ion-atom and ion-molecule reactions, 
the electron-density levels and distributions for the particular case reported 
here are not significantly affected. 

Reaction rate coefficients for reactions 14, 15 and 21 were taken 
from Ref. 34 . These rates are not well known but, fortunately, these re- 
actions do not have a significant influence on the results so that their un- 
certainty is unimportant compared with that of several other reactions in the 
model. 

In these calculations, the vibrational and electronic degrees of freedom 
are assumed to remain in thermal equilibrium with translation while the 
chemical reactions proceed at finite rates. The free- electron temperature 
was assumed to be equal to the heavy-particle translational temperature. 

On the basis of the in-flight measured electron temperatures presented in 
Section 5 of this paper, this assumption appears to be a reasonable one for 
the calculations reported here. 

2.4 Description of Experiment and Diagnostics 

The RAM C series of flights was designed to obtain data in the velocity 
regime of about 25, 000 ft/ sec where the importance of atomic particle 
ionization cannot be neglected. For the purposes of this paper, we are con- 
cerned with the second (RAM C-II) and final (RAM C-III) flights in this series. 
During this final flight, water and electrophilic liquid were periodically 
injected into the plasma layer to reduce the free- electron density level; how- 
ever, only data obtained in the absence of injection are presented. Many 
different diagnostic techniques were used on these vehicles, but we will 
place major emphasis on the electrostatic probe results here because this is 
the only diagnostic providing a direct measurement of electron-density profiles 
in the plasma. 

The reentry flight-path angle for these payloads was -15° and both of 
them were spin stabilized at 3 rpm. The payload geometries were nearly 
identical and consisted of a 6-inch radius hemispherical nose followed by a 
9° half-angle cone. The RAM C-II nose was covered by a beryllium-cap 
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heat sink during the portion of the e ^ tr ^ ablator for heat 

=s=sss;~ 

heat shield (Ref. 7). 

Figure 2 is a photograph of th. . C-U payload "gtte^locatiou 

of the electrostatic probe rake an a Vake so that the last probe extended 
were located on the leading edge of the rake so ‘hat the l P^ as 

to a maximum distance of approximately 7 cm. into the p lonaitudi „ al 

... . J TT'^ 71 -ru P rollectors were placed so that tneir longuuumat 

voltage so as to coueci ioxx». A t u_ data reduction procedure are 

construction, beyond that given m Fig. 3, and the ^ Teayxct p 
given in References 7. 11 and 12 and therefore will not be repeated n 

Figure 4 is a photograph of the RAM C-III payload illustrating the 

relative location of two rakes of electrostatic probes for tkls fll « ht ^ 
rake of constant bias-voltage probes was modified (Refs .2, 12) tkat ‘ h * 
number of collectors was increased to 16 and the distance that the final p 
extended into the plasma layer was increased from 7 to 14 cm. Other modi - 
cations of wedge and sweep angle were also made and are discussed m Ref . , 

12. The thermocouple rake was modified (Refs. 7,12) to include four vo ag - 
swept thin-wire probes located so that their longitudinal axis was approximately 
aligned with the flow direction. The outer wires extended to a maximum distanc 
of 9-5 cm. from the body into the plasma layer. The voltage sw ^F^ e * 
sweep cycle, and data collection procedure are all di sc uBSed in detatl m 

12 The relative locations of the probes and sketches of the rake con 
^ruction are given in Fig. 5 and 6 for the constant bias-voltage and swept- 
voltage probes, respectively. 

2. 5 Comparison of Theoretical and Measured Electron Densities 

The electrostatic probes discussed in the previous section were used 
to collect current, from which the electron density could be 

aUU^f^^^^'^o fhT bale^r Jit th^ probes continued 
to operatl and the voltages were recorded. However, for the purposes of this 
paper. ^ only Ihe electron-density data deduced for the altitude * 2 «. °°° 

to 275 000 ft will be compared with the theoretical resu . 
limitations of the theory, it was applied outside of what is considered the 
range of validity (due to the neglect of the *+/** term as a higher -order 
tJrm (Ref 15) in the present analyses) to calculate profiles for 214,000 ft. _ 
altitude Perhaps fortuitously, the calculated results compare favorably with 
the flight data athhis lower altitude as will be shown at the end of this section. 
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For the altitude range of interest here, the calculated mean free paths 
were such that the probes were operating in a free-molecular flow environ- 
ment. It was therefore possible to reduce the voltage- swept thin-wire probe 
data using the theory of Laframboise (Ref. 36) and the constant bias-voltage 
probe data using the theory of Smetana (Ref. 37). The accuracy of Laframboise' s 
theoretical results is well known. References 38-40 have demonstrated similarly 
good results using Smetana's theory. It can be shown that the number densities 
obtained using the results of Ref. 37 are insensitive to the magnitude of the 
electron temperature which is important from the viewpoint of flight-data 
analysis. This observation is not generally true when obtaining electron den- 
sities from probes aligned with the flow. 

Figure 7 illustrates the comparison of the calculated electron den- 
sities to those measured, on the C-II and C-III flight experiments, using 
the electrostatic probes previously described. Comparisons are made for 
altitudes of 233, 000 ft. , 250, 000 ft. , 265, 000 ft. , and 275, 000 ft. At all of 
these altitudes, the vehicle velocity was approximately constant at a value 
of approximately 25, 100 ft. / sec. The trajectories flown by these vehicles 
and their payload configurations were sufficiently similar that it was not 
necessary to perform separate calculations for the individual flights. There- 
fore, both the constant bias-voltage probe data of C-II and C-III and the 
swept- voltage probe data of C-III are included for comparison purposes. 

The uncertainty in the flight data obtained with the constant- voltage probes 
is indicated by the crosshatching and the vertical bars. The data reported 
are the time-averaged electron densities averaged over one body revolution 
and the indicated uncertainty represents the peak-to-peak fluctuations mainly 
due to angle of attack. A similar uncertainty existed in the swept- voltage 
probe data, but is not illustrated in order to prevent confusion. 

The results of the nonequilibrium flow, viscous shock- layer calcula- 
tions are indicated on Fig. 7 by the solid lines. The shock-layer thickness 
is also illustrated and shown to decrease with altitude up to approximately 
265, 000 ft. above which it doesn't appear to change much. The agreement 
between the calculated electron densities and the flight data is reasonably 
good over that portion of the shock layer for which data were obtained. The 
theoretical results show that the electron number-density levels increase 
near the outer edge of the viscous shock layer. Physically, this seems to be 
due to the diffusion of the chemical species into the shock-transition zone, 
giving a finite, non-zero level at the shock-layer edge. In addition, low 
temperatures exist near the shock-layer outer edge as a result of small shock 
angle at the base of the vehicle. This means a high level of mass density 
there, and, coupled with the finite electron species level, the electron number 
density is thus seen to increase near the outer-edge of the viscous shock layer. 
However, calculations were performed (Ref. 41) for a much longer cone (see 
Section 5 of this report) of the same nose radius and half-angle. The results 
show that the peak electron-density levels are located in the plasma layer away 
from the shock edge, partly due to the thinner, weaker shock existing there. 
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On the C-III flight the voltage- swept probes were used to make measure- 
ments of the electron temperature between 1. 5 and 9. 5 cm. from the body- 
surface. Since the translational temperature was used in the nonequilibrium 
calculations where, in fact, the electron temperature should have been used 
for evaluating many of the rate coefficient for reactions involving free electrons, 
these data are particularly important. The results of the electron temperature 
measurements are compared with the calculated heavy-particle translational 
temperature on Fig. 8. On the basis of this comparison, the electron tempera- 
ture appears to be approximately equal to the translational temperature at 
the measuring station. With the available flight data, it is difficult to be 
definitive about the relative magnitudes of the heavy-particle temperature and 
the electron temperature at other locations in the flow field. However, the 
flow fields of interest here are highly nonequilibrium and it can be seen from 
Table 1 that the reaction rate coefficients for the electron reactions are pro- 
portional to the -1. 5 to -4. 5 power of the electron temperature, suggesting 
that if there were a substantial difference between the inflight translational 
and electron temperatures, it would be reflected in the measured electron 
densities. It obviously cannot be stated with certainty that such is not the 
case. What can be said is that if one assumes these two temperatures to be 
in equilibrium throughout the flow field, then reasonably good agreement is 
obtained between data and theory at the measuring station. 

The calculations shown on Fig. 7 were obtained using the chemical 
model given in Table 1. The upper bound of the backward reaction rate co- 
efficients of reactions 7-10 and 13 were used in this calculation. The values 
used for the reaction rate coefficients for the remaining reactions were those 
given in the table. Because of uncertainties in rate coefficient values, it 
is important to assess their influence on the calculated electron densities 
and the subsequent agreement with flight data. This has been accomplished 
by repeating the nonequilibrium calculations using the lower bound of the 
backward reaction rate coefficients of reactions 7-10, and 13 while main- 
taining those for all other reactions at the Table 1 values. Figure 9 illustrates 
the results of these calculations for altitudes of 233, 000 ft. and 275, 000 ft. 

The influence at intermediate altitudes was similar to that shown. The un- 
certainty in reaction rate coefficients changes the predicted electron density 
by a factor of approximately two over the plasma layer. In general, the re- 
sulting uncertainty in the calculated number densities is consistent with the 
uncertainty in the flight data as can be seen from Fig. 9. 

It is also important to assess the influence of the chemical model on 
the calculated electron densities. Calculations were therefore performed 
for the probe location using one model in which the only positive ion was 
NO+ and a different model in which the positive ions were NO'*', N ? » 

N + , and 0 + . The results of these calculations are illustrated on Fig. 10 
for altitudes of 233, 000 ft. and 275,000 ft. The calculated electron densities 
at 233, 000 ft. were found to be essentially uninfluenced by ions other than 
NO+. However, at 275, 000 ft. inclusion of the additional ions in the chemical 
model drastically influenced the results. It is important to note that the re- 
sults of the calculation which included only the NO"* - ion fall well below the 
flight data in contrast to the more complete model which agrees well with 
flight data. At the probe location, the difference between the calculated results 
using these two models decreases monotonically in going from 275, 000 ft. to 
233, 000 ft. 


13 



The electrostatic-probe data were reduced assuming that the mass 
of dominant ion at the probe location was that of NO"^. Figures 11 and 12 
present calculated nonequilibrium species distributions across the plasma 
layer at the sphere-cone junction and at the probe location for altitudes of 
233, 000 ft. and 275, 000 ft. For the probe location and at 233, 000 ft. , the 
NO+ ion is dominant by approximately a factor of ten. However, at the 
junction the dominant ion was calculated to be 0 + . The calculations performed 
for the 275, 000 ft. trajectory point indicate that NO + was still dominant 
at the probe location but only by a factor of approximately two. The dominant 
ion at the junction was again 0 + but NO + and 0^ + were a close second. 

To test the limits of validity of the theory, calculations were performed 
for a trajectory point at 214,000 ft. At this altitude, conditions are such that 
theoretically the formulation should not be valid, since the Z^p/bX term was 
neglected in the present analysis as a higher-order term. However, the re- 
sults of the calculation are shown in Fig. 13 to be in reasonable agreement 
with the flight data. It is quite possible that this reasonably good agreement at 
214, 000 ft. may be fortuitous. However, the result may be of future interest 
if one is interested in obtaining an engineering approximation for similar 
vehicle and trajectory conditions. 

In addition to the electrostatic probe data discussed above, flush- 
mounted microwave reflectometers (Ref. 7) were also flown on the RAM vehicles. 
These particular reflectometers can be used to ascertain the electron-den- 
sity profile in the plasma layer under limited conditions, but they can generally 
be used to determine the peak value of the electron density at the antenna 
location. Figures 14 and 15 present favorable comparisons of the predicted 
electron- density profiles with the peak value of the electron density as deter- 
mined from the flight data (Ref. 9) for X/R N of 5 . 2 and 7 .65 (stations 3 and 4 
of Ref. 7) respectively, for altitudes of 233, 250, 265, and 275 kilofeet. Due 
to the nature of the reflectometers, plasma electron densities beyond ^/A 
of 0. 8 would not be seen. The agreement between the theoretical results 
and the flight data is considered to be reasonably good. 

Reflectometer measurements were also performed in the nose region 
at X/ R n of about 0. 9 and just after the hemisphere-cone junction at 2. 1. 
However, the agreement between theoretical prediction and flight data was 
not nearly as good for these locations as it was for the two downstream 
stations. In the nose region, the calculated peak number density was about 
20 times greater than the reflectometer data and at X / R.^ of 2. 1 the calcu- 
lated peak value was about 10 times greater than the reflectometer data. 

At the present time, the reason for the discrepancy between the 
reflectometer data at X/R^ of 0.9 and 2.1 and the theoretical prediction is 
not understood. As possible causes, the translational- vibrational temperature 
nonequilibrium and the translational-electron temperature nonequilibrium were 
investigated for an altitude of 265,000 ft. A particular streamline was chosen 
which enters the viscous shock layer at about 30 degrees from the stagnation 
line. The results of the calculations of the vibrational relaxation distance for 
molecular nitrogen indicated that translational -vibrational equilibrium was 
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achieved in less than 0.8 cm. behind the shock. Molecular oxygen would 
achieve translational-vibrational equilibrium in significantly less distance. 

It therefore appears reasonable to assume translational -vibrational equilibrium 
for the viscous hypersonic flow around a RAM-C vehicle -at high altitudes. 

In addition, comparison of the distribution of the electron temperature along 
the streamline noted above with the nitrogen vibrational temperature and the 
heavy-particle translational temperature showed that it is appropriate to 
assume translational-electron temperature equilibrium in the nose region. 


2. 6 Conclusions 

In this paper, a multicomponent, ionized, viscous shock-layer flow 
about a blunt-nosed vehicle (sphere- cone) at high altitudes has been theo- 
retically analyzed. This was accomplished by transforming the coordinates 
and the flow variables, and by utilizing the Karman-Pohlhausen integral 
method to the Navier-Stokes equations and the species-conservation equations 
under the thin shock-layer assumption. 

The theoretical development was applied to a RAM-C sphere-cone 
vehicle (9° semi-cone angle) developed by NASA because of the wealth of 
in-flight measurements. Eleven chemical species (CR, N , O, N, NO, 
NO+, 0+, N + , 0 2 + , N~ + and e“ ) were included in tne airalysis for eleven 

reactions involving only the neutral species and fifteen reactions for charged 
species. Solutions were obtained at various altitudes where the viscous 
shock-layer flow exists around the vehicle. The theoretical results for the 
electron-concentration profiles compare favorably with the experimental 
data obtained using electrostatic probes located at the base of the vehicle. 

In addition, the heavy-particle translational temperatures obtained from 
theory seems to agree fairly well with the electron temperatures measured 
from the voltage swept thin-wire electrostatic probes. 

The effect of the reaction- rate coefficients on the quantitative theo- 
retical results was also analyzed and shown to be within the uncertainty range 
of the flight data. The results also demonstrate the importance of including 
in the chemical model the positive ions N'f, O^", ac ^ : *-tion 

NO+ for high altitudes and velocities considered here. 
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3. ANTENNA IMPEDANCE DETERMINATION OF ELECTRON DENSITY 


3 . 1 Introduction 


Many (Refs. 7, 9 i 42-48) different diagnostic techniques have been em- 
ployed to determine the electron-density level in the plasma layer surrounding 
an entry vehicle. Several (Refs. 9 >42-48) of these instruments can be located 
flush with the vehicle skin while others (Ref. 7) protrude into the flow field. It 
is generally recognized that if one is interested in obtaining information about 
the electron-density gradient and level in the plasma layer, then the protruding 
probes are the most desirable. However, for some vehicles and under certain 
flight conditions, protruding objects are not feasible thus requiring use of flush- 
mounted diagnostics in order to obtain as much of the desired information as 
pos sible . 

Bailey and Swift (Ref. 49) have treated theoretically the problem of a 
circular - ape rture antenna excited in the TE j j mode. Swift, et al . (Ref. 46) and 
Grantham (Ref. 9) have used circular antennas on a re-entry vehicle to obtain in- 
flight admittance measurements from which they could deduce the plasma prop- 
erties in the near field. 

The accuracy of the calculational technique (Ref. 49) has not been assessed 
for a flowing plasma and it is therefore the purpose of this study to investigate 
the agreement between experimental and theoretical antenna admittance for a 
known distribution of electron density in the plasma layer adjacent to a flat- 
plate surface located in a flowing plasma. The experimental phase of this work 
is an extension to the previous work of Taylor (Ref. 50) who utilized the same 
antenna employed here to perform admittance measurements in a shock-tube 
flow but under much different gas-dynamic conditions. His boundary layers 
were significantly thinner than the present ones and the axial length of the 
uniform plasma was much less. The collision frequency, of his plasma 

(0.1 torr experiments) was less than 0.15 which is of the same order as that 
obtained here. Taylor illustrates good agreement between his admittance data 
and the predictions of Croswell, et al,(Ref. 5 1). However, because his boundary 
layers were thin, it was not possible to probe them in order to obtain a 
measured electron-density profile for use in the calculations. 

3. 2 The Input- Admittance Method of Plasma Diagnostics 

The theory of the input admittance of plasma-covered aperture anten- 
nas has been studied by a number of investigators over the past several years. 
These studies were motivated by the desire to extend the diagnostic range of 
antennas beyond a simple marker of critical electron density. If the antenna 
consists of a waveguide opening onto a flat ground plane, the normalized 
input- admittance expression is given by 
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where V c is the characteristic admittance of the waveguide, 
impressed voltage, is the electric field at the aperture, 

netic field, & S is the surface element in the plane, and 
the ground plane. The relationship between input admittance 
reflection coefficient T is: 


r = 


i -x, 
1 + X* 


_V 0 is the 
H 5 * is the mag- 
Z = 0 defines 
and the complex 


( 2 ) 


The near-field properties of the circular- wave guide antenna were 
investigated by Bailey and Swift (Ref. 49)* and later generalized to account for 
stratified layers of inhomogeneous dielectric material (including plasmas). 

The formulation of the problem resulted in a reduction of Equation (1) to a 
single integral which was numerically evaluated as a function of collision 
frequency, electron density, and gradients of electron density. 


Prior to the shock-tunnel experiments, the input admittance was cal- 
culated as a function of boundary-layer thickness ( 6 ) and free- stream 
electron density for the antenna shown in Fig. 16. The plasma profiles used 
for these calculations (given in Fig. 17) were assumed to be parabolic, with 
parameters adjusted so that the electron density was zero at the wall and 
the free- stream value at 8 . The particular values of 8 selected for the pre- 
experiment analysis were 0, 0.2, 0.4, 0.6, 1.2, and 2.4 inches. Beyond 
the point at which n c /n e ) cr reached 1, the plasma was assumed to be homoge- 
neous to 10 cm which was sufficiently thick so that reflections from the air- 
plasma interface at 10 cm were negligible. 

The theoretical results are shown in Fig. 18 as a function of <f and 
n e./ rl t ) cr for V/tt) = 0. 1. These results suggest the following diagnostic 
applications cf the circular-aperture antenna: 

1. The measurements of the initial response of the antenna near 
r 'e./ n e' l c .r = 1 can be used to uniquely infer both n e /n e ) cr and 8 . 


2. When the plasma is overdense, the admittance measurements 
can be used to infer either n ft /n e ) cr or € . Some "a priori" knowledge of 
one of the parameters must be given in order to provide a measurement of 
the other. 

The boundary-layer electron-density profile was measured at two 
locations on the flat -pi ate model for one of the experimental conditions used 
in this work. With this information, it was possible to theoretically estimate 
the magnitude of the boundary-layer thickness for the remaining experimental 
conditions. By assuming that the peak electron density occurs at the boundary- 
layer edge, which introduces some error because previous experimental 
results (Refs. 52,53) for similar conditions suggest that it occurs at distances 
slightly less than the boundary-layer thickness, it was possible to obtain reason- 
able estimates of electron density in the overdense region. It is worth noting 
that the antennas and the calculations reported in Fig. 18 can be scaled for 
application at frequencies other than S-band. The antenna for which these 
calculations and experiments were performed is a reproduction of one flown on 
a re-entry experiment (Ref. 8) and will henceforth be designated as the RAM C- 
III S-band antenna. 
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Another type of circular-aperture antenna was also used in these 
laboratory experiments and in the RAM flight experiments. This antenna 
was designed for and used on the RAM C-II flight two years before the 
boundary- value problem of the circular- waveguide antenna was developed. 
Although the antenna operated at the same frequency as the one discussed 
above, there were some major design differences. Referring back to Fig. 16 » 
these differences are: 

1. The antenna was extended through the dielectric material so 
that the aperture was contiguous to the plasma. 

2. The diameter of the antenna was designed to be smaller than 
that of the one shown in Fig. 16 in order to reduce the effects of lateral gra- 
dients of electron density over the aperture. 

3. Because the waveguide section was near cut-off, tuning screws 
were inserted to match the antenna to free space. This tuning was accounted 
for in the theoretical model. 

Subsequent to the flight of RAM C-II, the aperture admittance was 
computed using Equation (1) in connection with the transformation derived in 
Appendix A. This transformation was used to calculate the admittance in 
preparation for the tunnel experiments described herein. 


3. 3 Experimental Apparatus and Technique 

An engineering drawing of the RAM C-III S-band antenna used in 
the experiments is shown in Fig. 16. It consists of a circular cavity 4. 5 
inches deep and 2. 65 inches in diameter. A co-ax-to- waveguide transition 
was inserted 0. 4 inches from the back wall of the cavity, i. e. , one-quarter 
of a guide wavelength. The length and diameter of the probe were adjusted 
to obtain a 1:1 impedance match between the waveguide and RG-141 coaxial 
cable at the operating frequency of 3. 348 GHz. An 0. 8 inch teflon plug was 
inserted into the waveguide flush with the aperture and the ground plane was 
covered with 0. 3 inches of plexiglas. The reasons for selecting these dimen- 
sions is discussed elsewhere (Ref. 46). 

The antenna described above and a second one mentioned previously 
were mounted in a sharp flat-plate model which was subsequently placed in 
the fiberglas nozzle of the shock tunnel (Ref. 52). The C-III antenna was operated 
at 6. 75 inches from the plate leading edge and the C-II antenna was operated 
at 11. 25 inches from the leading edge. The boundary layer over the plate 
provided an electron-density gradient in the near field of the antenna. This 
b oundary- laye r was sufficiently thick so as to permit good resolution of the 
electron-density profile by detailed probing with thin-wire Langmuir probes. 

The flat-plate model was 7-inches wide by 30-inches long. The top 
surface was constructed of plexiglas with the exception of a 1-inch long stain- 
less-steel leading edge. Boundary- layer surveys were performed at distances 
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from the leading edge of 6. 75 and 11. 25 inches. At the most upstream location, 
the flow was considered to be two-dimensional since spanwise static-pressure 
measurements obtained by Boyer, et al . (Ref. 53) at 5 and- 10 inches from the 
leading edge, for similar flow conditions, indicated that the flow was uniform 
across the plate. However, for the 11.25 inch location pressure measurements 
were not obtained, and it is possible that the flow here may hot have been two- 
dimensional. 

During all of the experiments noted above, an independent measure- 
ment of the integrated free-stream electron density was obtained at 11 inches 
and one inch upstream of the flat-plate leading edge using microwave inter- 
ferometers operating at frequencies of 35 and 17 GHz, respectively. 

J 

The probes used in these experiments were constructed by surrounding 
0. 004 inch diameter tungsten wires with a quartz envelope, leaving a nominal 
0. 400 inch length of bare wire exposed. Immediately prior to each run, the 
tungsten oxide was removed by placing the probe in a dilute solution of sodium 
hydroxide and passing approximately 400 (xa of current through the circuit 
for approximately 10 minutes. In performing the probe measurements, care 
was taken to keep the probe holder far out in the flow to avoid any possible 
interference effects in the subsonic-flow region. The quartz tubing surrounding 
the probes was bent as shown in Ref. 5.2. The 1/16-inch diameter tubing used 
in the probe construction did not permit probing to distances less than 0. 060 
inches off the plate surface. 

Figure 19 illustrates the two S-band antennas mounted in the flat-plate 
model. Also shown is an X-band horn antenna, the results of which are not 
reported here. A schematic diagram of the four-probe reflectometer is given 
in Fig. 20. The reflected signal was also monitored as indicated on the sche- 
matic. This system was convenient for data reduction purposes in that five 
arcs were always available for Smith chart construction. Waveguide instead 
of coaxial cable was used wherever possible in the system. 

As part of the pre-experiment calibration, the diode spacings were 
measured using a moveable short. If the microwave diodes are numbered so 
that #1 is nearest the generator and #4 is nearest the short, then the spacings 
were 1-2 = 0.-1255 A,., 2-3 = 0. 1218 A c , 3-4 = 0. 1278 A fl , and 4-1 = 0. 1252 

where was the guide wavelength of 11.378 cm. The reflectometer system 

was tuned using a matched load with VSWR of approximately 1. 025 to 1. 035. 

Eachof the diodes was then calibrated as follows. An attenuation of 3 db was 
set on the variable attenuator and the four diode outputs were adjusted until 
they were equal. The attenuation was increased in increments of 1 db until reaching 
26 db. The attenuator was then set to o © attenuation and the final readings 
taken. The voltages were normalized so that the diode output at 3 db attenua- 
tion represented the diameter of the Smith chart. During the experiments 
this 3 db tare was always set on the variable attenuator. 

The free- space admittance was measured with the antenna inside and 
outside of the shock tunnel, but no difference was observed between the re- 
sults. Aluminum foil was then placed over the aperture of the antennas 
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obtain the antenna- short reference. In addition, the coaxial cable was shorted 
at the antenna input terminal to obtain an additional reference. 

The antenna admittance and the Langmuir -probe measurements were 
performed in separate experiments. In both cases, the test gas was carbon 
monoxide but the level" of residual nitrogen present in the test g&s was dif- 
ferent. The carbon monoxide used in the antenna experiments was supplied 
by Lif-O-Gen, Inc. A chemical analysis of the gas indicated the following: 

25 ± 5 ppm of nitrogen, argon less than 2 ppm, carbon dioxide less than 

10 ppm, helium less than 2 ppm, oxygen less than 2 ppm and water less than 
1 ppm. Unfortunately, a sufficient quantity of this gas could not be obtained 
for the boundary-layer measurements and it was necessary to use ultra 
high-purity grade supplied by Air Products and Chemicals, Inc. The prin- 
cipal impurity in this gas was nitrogen estimated to be on the order of 1500 
to 2000 ppm. The presence of this additional nitrogen was important from 
a chemical kinetics viewpoint and had to be considered. However, even with 
the nitrogen, the dominant ion in the free- stream flow was calculated to be 
C + consistent with the calculations for the purer gas. 


3.4 Electron-Density and Electron-Temperature Measurements 

The boundary- layer measurements of electron density and electron 
temperature reported here were performed in carbon monoxide that had 
expanded from equilibrium reservoir conditions of 7060° K at 17. 3 atm. 

The relative magnitudes of the mean free paths A;. n , ^n-n ? A e -n v Kx.-* , 
and A e _ e were calculated using the expressions summarized by Sonin (Ref. 5. 
and presented in Fig. 18. The electron depsity was obtained from the ion-current 
portion of the Langmuir-probe curr ent-voltage characteristic using the free- 
molecular flow theory of Laframboise (Ref. 36). In some cases the electron- 
current portion of the characteristic could be used to obtain electron density which 
was found (Ref. 55) to be in good agreement with the value obtained from the ion 
current. An ion mass of 12 was used in the data reduction consistent with the 
calculated species concentrations presented in Ref. 30. 

Electron-density measurements obtained in the boundary layer at 
6. 7 5 inches from the plate leading edge are shown in Fig. 22. The measured 
boundary- layer values approach the estimated free- stream value at approxi- 
mately 1. 2 inches from the surface. This estimate of the free- stream number 
density was obtained from microwave-interferometer data taken at the measuring 
station in the absence of the plate, assuming that the plate does not significantly 
disturb the flow. 

Since the dominant ion in the free- stream flow was C + and the electron- 
temperature measurements shown in Fig. 23 suggest a relatively constant 
value, it is unlikely that gas -phase recombination of ions and electrons would 
be important. However, it is difficult to be certain as to whether or not the 
plexiglas surface would provide a good catalyst for surface recombination 
of ions and electrons. If the wall were catalytic, then a diffusion controlled 
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boundary layer would be appropriate and it would have a number -density 
profile consistent with the measurements. One cannot be certain without 
detailed calculations which were beyond the scope of this work, but it is 
likely that the boundary-layer gas phase chemistry was frozen at an upstream 
value and that surface recombination controlled the boundary- layer electron- 
density profile. 

The boundary- layer electron-density and electron-temperature pro- 
files obtained at 11.25 inches from the leading edge are presented in Figs. 22 
and 23. The number-density measurements suggest a slightly thicker boundary 
layer at this location but the profile shape is similar to that measured at 
6. 75 inches. The electron temperature was again found to remain rather 
uniform through the boundary layer and along the plate. 

The antenna measurements to be discussed later were performed for 
several reflected- shock reservoir conditions in order to achieve a range in 
the magnitude of ^e/^ e ) cr » but the boundary- layer measurements were per- 
formed for only one experimental condition because of the detailed nature of 
this experiment. The range of test conditions was such that boundary- layer 
electron-density profiles in terms of ij/rie),* vs ij_/S could be constructed 
and assumed similar. This profile could then be appropriately adjusted for 
the different free-stream conditions. To aid in constructing the electron- 
density profiles for conditions other than the measured one, a simple estimate 
of the boundary- layer thickness, 6 , was required. The technique selected 
was that described by Burke (Ref. 56) and Figs. 22 and 24 show that the results 
of this method are somewhat less than the value that one might have estimated 
on the basis of the electron-density measurements . However, for the intended 
purpose of scaling the profiles, this result is quite satisfactory. 


3. 5 Admittance Measurements 

The RAM C-III antenna was installed in the flat plate at 6. 75 inches 
from the leading edge and the top side of the plate was located on the nozzle 
centerline with the leading edge at 22. 5 inches from the nozzle throat. For 
these experiments, the peak electron density over the antenna was varied 
in the range of approximately 0. 5 < / n e) cr 5 by adjusting the strength 
of the incident shock. The experimental results obtained, including the arc 
lengths for the five diodes, are summarized in Table II. 

The RAM C-II antenna was operated at 11. 25 inches from the leading 
edge. The major differences between the 6. 75 and 11. 25-inch measuring 
stations are the boundary- layer thickness and the peak value of the electron 
density. Several experiments were performed at the 11,25-inch location 
and these spanned the range from 0. 2 < / Kl e ) cr ^ 2 . 1 . The results of these 

experiments are given in Table III. 

Figure 26 is typical of the oscilloscope records obtained for the reflec- 
tometer diodes. The test-flow duration is noted on the oscilloscope records. 
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The 50-u.sec duration pulse near the right side of these records is the zero- 
signal level obtained by turning off the TWT power supply with a gating 
pulse. The horizontal line is the no-plasma signal level recorded approximateb 
two milliseconds prior to plasma arrival at the plate leading edge. A limited 
number of lower confidence -level data points can be obtained from the post 
test-time plasma. In some cases this was done and the data points are so 
designated in Table II and III. 

All of the reflectometer data obtained in the CO plasma with this 
antenna are compared with the theoretical predictions on Fig. 27 for the 
RAM C-III antenna and on Fig. 28 for the RAM C-II antenna. The circles 
represent the uncertainty of the measurements as determined by the inter- 
section of the five diode voltages. Much of the phase rotation uncertainty asso- 
ciated with these records was due to the #1 diode arc being consistently long. 
Care was taken to calibrate each of the oscilloscopes just prior to and after 
each experiment using a vacuum-tube volt meter, thus it is doubtful that the 
difficulty is associated with a calibration error. This diode appeared to 
operate satisfactorily when obtaining the free- space admittance for which the 
circle of uncertainty was comparatively small. The estimated uncertainty 
in the value of "e. / is also g iven on Figs. 27 and 28. 

The theoretical results presented in Fig. 27 were obtained using a 
parabolic profile of electron density that represented the best fit to the 
el ect ron -density profile measured at the antenna location (8^~'1.2 inch). The 
comparison between the theoretical and the experimental values of n g/ n e^cr 
presented in Fig. 27 indicates a consistent difference of about 20 to 40 percent. 
This discrepancy is possibly due to an uncertainty in the experimental deter- 
mination of the free-stream electron density over the antenna. Recall that 
the admittance measurements and the boundary-layer probing measurements 
were performed in separate experiments in order to avoid interference effects. 
Microwave interferometer measurements of the averaged free-stream electron 
density were performed at 1 inch and 11 inches upstream of the plate leading 
edge for all experiments and thus provided the only quantitative link between 
the measurements. When obtaining the admittance measurements, these 
microwave interferometer results were extrapolated to the antenna locations 
in order to obtain an estimate of the free-stream electron density assuming 
that the expanding plasma was not influenced by the presence of the model. 
Figures 22 and 24 illustrate that the extrapolation described above is reason- 
ably good, but an error of 20 to 40 percent in the estimated electron density 
would be within the uncertainty. 

Figure 28 presents a comparison of the experimental data with the 
theoretical results for the RAM C-II antenna. These results were obtained 
with the antenna located at 11.25 inches from the plate leading edge. The 
experimentally determined electron-density profile given in Fig. 24 was used 
in order to obtain the theoretical results presented in Fig. 28. The results 
obtained using this antenna are similar to those obtained for the C-III 
antenna in that the experimentally determined reflection coefficients are about 
30 percent less than the theoretical values at values of n e /n e ) cr less than 1. 
The agreement becomes better as the electron density ratio increases to 
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values greater than 1. The experimentally determined phase rotations are 
shown to be generally in better agreement with the theory except at the very 
small values of n e /n e ) cr . Phase measurements at low values of n e /n e ) cr 
is more difficult than at higher values because of the small reflection 
coefficients . 


3. 6 Conclusions 

Two circular-aperture flush-mounted antennas, for which a theo- 
retical description of the antenna- plasma interaction has been developed, 
were operated in the presence of ionized boundary layers of low collision 
frequency. The boundary- layer electron-density profile was independently 
measured using voltage-swept thin-wire probes and this measured profile 
was used as an input condition for the theory. Both antennas demonstrated 
marked increases in reflection coefficient and phase rotation when the peak 
electron density over the antenna approached critical. In general, good 
agreement was obtained between the predicted and measured reflection coef- 
ficient and phase rotation over the range of /n e ) er studied. These results 
illustrate that the input admittance method can be effectively used to obtain 
local measurements of electron density. Calculations of the admittance of 
a flush-mounted antenna have been presented as useful working curves from 
which boundary- layer thickness and free-stream electron density can be 
determined. These curves can be used at arbitrary frequencies by appro- 
priate electromagnetic scaling. 
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4. USE OF POSITIVELY BIASED ELECTROSTATIC PROBES TO OBTAIN 
ELECTRON DENSITY IN COLLISIONLESS FLOWS 


4,1 Introduction 

There are several well known (Refs* 57-61) advantages of using the 
electron- cur rent portion of an electrostatic-probe cur rent-voltage character- 
istic for determination of electron density. Perhaps the most important 
advantages are that the collected current is independent of the ion mass and 
also independent of the probe orientation relative to the mass motion of heavy 
particles because of the relatively high electron thermal speed. In addition, 
the collision cross sections for electron-particle interactions are generally 
smaller than those for ion-particle interactions suggesting that the collisionless- 
flow approximation should be valid over a greater range of gas densities. 
Further, the electron current collected by the probe can be made insensitive 
to the electron temperature by working at large probe voltages. Sutton (Ref. 62) 
has described the feasibility of using thin-wire probes biased for electron 
current in order to measure electron-density fluctuations in hypersonic turbu- 
lent wakes. However, relatively few (Ref. 61) experimental studies have been 
reported in the literature that have attempted to show that the electron-current 
region can be successfully used to obtain these electron-density measurements 
within the framework of existing theoretical results. It is thus the purpose of 
this paper to report the results of an experimental study undertaken in order 
to determine if the existing theoretical results of Laframboise (Ref. 36) for the 
current collected by thin-wire probes operating in the free -molecular flow 
regime could be successfully used in the electron-current region. Many people 
have illustrated the correctness of this theory for the ion-current region. The 
details of Laframboise' s work are well known and will not be described here. 

The current studies were motivated by the availability of in-flight 
voltage - swept probe data obtained by Jones (Ref. 63) as part of the Langley 
Research Center RAM program. The voltage applied to his probes was 
sufficiently positive that the probes were swept well into the electron-current 
region at altitudes for which the free -molecular flow approximation should be 
valid. Electron densities have been obtained from the ion-current portion of 
his probe characteristics and are reported elsewhere (Ref. 63). 

Center (Ref. 61) used a pressure-driven shock tube to produce a plasma 
behind the incident shock in which he performed probe experiments under free- 
molecular flow conditions. The voltage on his probes was held constant during 
the experiment because of the rapidly varying conditions behind the incident 
shock. By performing a series of experiments he was able to construct the 
current-voltage characteristic and subsequently show that the experimental 
data are predicted by the theoretical results of Laframboise (Ref. 36) . The 
results presented in this paper are in agreement with those described by 
Center. However, the fluid-flow environments in which the experiments were 
performed were significantly different and the probe characteristics were 
obtained by sweeping the probe voltage instead of maintaining constant voltage. 
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4.2 Experimental Apparatus and Technique 

A pressure-driven shock tube was used to produce a reservoir of 
high-temperature gas which was subsequently expanded in a conical nozzle 
constructed of Fiberglas. A detailed discussion of this shock tube and nozzle 
is given in Refs. 64 and 65. 

The majority of the probe measurements reported here were performed 
in the boundary layer of a sharp flat plate mounted in the nozzle flow such that 
its leading edge was 22. 5 inches from the throat and its top surface was on 
the nozzle centerline. A few experiments were also performed in the absence 
of the plate with the probes located on the nozzle centerline at 32. 5 inches 
from the throat. An independent and simultaneous measurement of averaged 
free-stream electron density was obtained at 11.5 and 21. 5 inches (31. 5 inches 
for the free-stream experiments) from the throat using microwave interfero- 
meters operating at 35 and 17 GHz, respectively. Carbon monoxide was used 
as the test gas for the boundary- layer experiments. The carbon monoxide 
test gas was expanded from an equlibrium reservoir condition of 7060° K at 17.3 
atm pressure. However, the free-stream measurements were performed in 
both carbon monoxide and nitrogen test gases. The nitrogen test gas was 
expanded from an equilibrium reservoir condition of 7200° K at 17. 1 atm pres- 
sure. Previous boundary -layer experiments have been reported (Ref. 52) which use 
nitrogen as the test gas but in these measurements the probes were not swept 
into the electron-current region. 

The carbon monoxide used here was ultra high-purity grade supplied 
by Air Products and Chemicals Inc. The principal impurity in this gas was 
nitrogen estimated to be on the order of 1500 to 2000 ppm. The nitrogen test 
gas was ultra-pure carrier grade also supplied by Air Products and Chem- 
icals Inc. A chemical analysis of the nitrogen gas indicated the following: 
oxygen less than 0. 5 ppm, total hydrocarbons less than 1 ppm, and water 
less than 0. 15 ppm. 

The probes used in these experiments were constructed by surrounding 
0. 004 in. diameter tungsten wires with a quartz envelope, leaving a nominal 
0. 400 in. length of bare wire exposed. Immediately prior to each run, the 
tungsten oxide is removed by placing the probe in a dilute solution of sodium 
hydroxide and passing approximately 400 pa of current through the circuit 
for approximately 10 minutes. 

4. 3 Discussion of Results 

By using Laframboise 1 s (Ref. 36) theory to obtain electron temperatures and 
electron densities from the experimentally determined probe characteristics, 
it has been assumed that the flow conditions were collisionless with respect 
to the probe diameter. At all of the measuring stations, the electron mean 
free paths were many times greater than the probe diameter. Because of the 
relative insensitivity of electron collection to the magnitude of the ion-neutral 
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and neutral-neutral mean free paths, it is anticipated that the electron tempera- 
ture deduced from the electron- retarding region and the electron density- 
deduced from the electron- current region of the probe characteristic should 
be correct. The ion- neutral mean free path was the shortest, being slightly 
greater than the probe diameter at the upstream measuring station. How- 
ever, the electron-density results presented later in this paper suggest 
that the ion- neutral mean free path apparently had a relatively insignificant 
influence on the collected ion current. The calculated mean free paths and 
species concentrations in the expanding carbon monoxide and nitrogen flows 
used for the present experimental test conditions are given in detail in Refs. 30 
aid 66, respectively. 

It is important to note that the carbon monoxide test gas used in the 
work reported in Ref. 30 had substantially less impurities (25 ppm nitrogen) 
than that used in the present experiments (1500 to 2000 ppm nitrogen). It 
was not possible to obtain additional CO from Lif-O-Gen Inc. with the low 
nitrogen content. The additional nitrogen appeared to influence the chemical 
kinetics of the expanding plasma by decreasing the electron-density level by 
a factor of approximately two throughout the expansion. For the purposes 
of this study the details of the kinetics are important and had to be considered. 
Even with the nitrogen present, the dominant ion in the expansion was calculated 
to be C + so that ion mass 12 was used inlhe analysis of the ion-current data. 
The method of probe-data reduction used here has been outlined in detail in 
Ref. 67 and will not be repeated. The familiar normalized ion-and electron- 
current densities( .+ + and j, _ ) were taken from Laframboise's theoretical 
results for / T e = 0 and are given in Figs. 29 and 30 .respectively. 

For the experimental conditions considered here the Tx,/T e - 0 theoretical 
results are satisfactory since the free-stream ion temperature to electron 
temperature ratio was on the order of 0. 1 to 0. 2. 

Table IV presents a summary of the experimental results obtained as 
part of this study. Included in this table are the electron densities deter- 
mined from the ion- and el ectron- cur rent portions of the current- voltage 
characteristics, the calculated values of R^>/ A D for these regimes, a 

comparison of (.V^o _ ) obtained from the tangent intercept technique and 

an estimate obtained from the theory, the electron temperatures determined 
from the electron- retarding region, and electron densities determined using 
a microwave interferometer when free-stream probe measurements are given. 

Figures 31 and 32 present the ion-current region and the electron-re- 
tarding and electron-current regions respectively of the probe characteristic 
obtained in the boundary layer for a typical experiment. On Figure 31 the ion 
current is utilized to deduce the number-density at ”2.5, For this par- 

ticular case the value of (R^ / A p ) 2 J-+ was 16. 55 which is in the orbital- 
motion limited regime. The corresponding value of R^/A 0 was l* 72 and the 
floating potential was -0. 43 volts. The electron density calculated using 
the theoretical results of Laframboise (Ref. 36) was 1.2 x 10 e /cm . Also 
shown on Fig. 31 is the electron density calculated using the electron current at 
potential of % io = +5. The electron current at % ^ =+7.5 was also used but 
the resulting electron density was essentially unchanged. The value of 
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(R p /A D ) X. at X ^ = +5 was calculated to be 7. 36 and the corresponding value 
of was 1.78 which is reasonably close to the value calculated from 

the ion current. The number density deduced from the electron current was 
1.3 x 10l0 e”/cm^ which is in very good agreement with the value previously 
obtained from the ion-current portion of the characteristic. 

The electron temperature was obtained from the results presented in 
Fig. 32. The value of e. /&T e for this particular experiment was 5. 37 or 

T e = 2160° K. By extrapolating the line faired through the retarding-field 
region to intersect with a line faired through the electron-current region 
one can obtain an estimate of the plasma potential. For this case, the extrapo- 
lation gave a plasma potential of +0. 5 volts; the floating potential, which is 
easily determined from Fig. 31 was -0.43 volts resulting in an experimental 
value of ( - Vf ) =0. 93 volts. Using the theoretical expressions given 

in Eqs. 1 and 2 of Ref. 59 and equating to at the floating potential, one 
can obtain an estimate of the value of (V» - from the theory to be given 

by ( /e ) Ln ( rr\^ / ) 0,5 j. + . We will refer to this value as ( 

and compare it in Table IV with (V* - Vp) ei6> determine from the experimental 
data using the tangent intercept method in a manner illustrated on Fig. 32. 

For the purposes of this paper, was assumed to be equal to 1. 0 which 

results in a maximum error in (V^ - Vf ) of 10 to 12 percent for the conditions 
of these experiments. 

Assuming C + to be the dominant ion, then (V^ - V f ) est — ( T e / 2320 ); 
then assuming N + to be the dominant ion in the nitrogen experiments^ gives 
(\4»~Vp) e st S T e /2290. Using the previous expression for the C + ion 
the value of (V„ - Vf ) e -st for the results presented in Figs. 31 and 32 is equal 
to 0. 932 as compared to the experimental value of 0. 93. The agreement 
between (V^- estimated from the theory and the value obtained from the 
experiment was not always this good (as illustrated in Table IV but in general 
the agreement was within 0.2 volts. 

Table IV also contains number-density results obtained in the free- 
stream flow using both the microwave interferometer and voltage- swept probes. 
For the CO plasma, the probe-determined electron densities were in good 
agreement with each other (ion and electron current) but in two of the three 
cases they were about 30 to 40 percent less than the microwave results. This 
same general trend was reported in Ref, 30 but the reason for the disagree- 
ment is unknown. Probe and microwave-interferometer results for a nitrogen 
plasma are also reported in Table IV. With the exception of the final entry in 
the table, the number densities determined from the ion- and electron-current 
portions of the characteristic are in good agreement. In all cases, the number 
densities determined from the ion current were in good agreement with the 
mic rowave- interferometer measurements. 

Figures 33 and 34 present a comparison of the experimentally determined 
cur rent- voltage characteristic with the theoretical characteristic of Laframboise. 
The data used here are the same as those given in Figs. 31 and 32. To obtain 
the comparison, the plasma potential ( ) was assumed to be 0. 5 volts con- 

sistent with Fig. 32 and the experimental characteristic was matched to the 
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have been made wholly or in part in the previous analyses and the justification 
for these assumptions can be found in Refs. 15, 16, 18, 19. 


Based on the above assumptions the governing equations for the viscous, 
ionized shock layer in the merged-layer regime become 





t ‘ (6) 

These equations express the conservation of mass, momentum, energy, and 
chemical species, respectively, along with an equation of state* (See Fig* 1 
for a description of the flow field* ) The eleven chemical species considered 
are: C> 2 , N 2> O, N, NO, NO+, O+, N+, C> 2 + , N + and e“ . With 
the two-layer model of Cheng, Equations 1-6 are appued to the viscous 
shock layer allowing for the diffusion of the species into the shock-transition 
zone by using the modified Rankine-Hugoniot conditions. Analysis of the 
entire flow field as performed in Ref. 19 reinforces the applicability of the 
two-layer model, and thus only the viscous shock layer will be considered 
here. The chemical reactions and their associated reaction rate coefficients 
used simultaneously for dissociation and ionization are given in Table I and 
are discussed in detail in Ref. 69* It is appropriate to note here that eleven 
reactions are used involving only neutral species and fifteen reactions in- 
volving charged species. 
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In obtaining solutions to the differential equations, the coordinates 
and the flow variables were first transformed, and the Karman-Pohlhausen 
integral method was applied because of the ease of application and the rela- 
tively sm 9 .ll computation time without unduly sacrificing the accuracy o e 
results (Refs. 2,21,69)* The transformed equations were integrated along the 
body from the stagnation region through the sphere -cone junction to the afterbody 
conical section by using the Adams-Molton predictor-corrector method. 

In order to simplify the analysis while accounting for changes in shock shape 
along the body, the distribution of the shock angle was initially assumed and 
updated from the results. Based on the cases considered in the present 
analysis, two iterations were usually sufficient to obtain the correct s 
shape for the body. The solutions are obtained in terms of the flow variables 
such as the velocity and temperature, and in terms of the chemical- species 
distributions in the shock layer and along the body surface. 


The present analysis was applied to two different bodie s: a 9° sphere 
cone with a 6 inch nose radius at 233 kft and a 20° sphere cone with a 4 ft nose 
radius at 280 kft and 310 kft. The latter geometry was used to approximate 
the windward plane of symmetry flow configuration of a space _ shuttle - | ® 
free-stream velocity for both bodies was taken to be 25, 100 ft/ sec. Solutions 
were obtained for a flow field including not only the nose region but also the 
far downstream region (.%/ R N 90) . 


5. 3 Results and Discussion 


Figure 35 shows the distributions of the viscous shock-layer thickness 
and the temperature behind the shock for the two (9° and 20°) sphere-cone 
bodies. It is noted that the temperatures behind the shock m the nose region 
are about equal for both bodies, due to the same free- stream kinetic energy 
being converted into thermal energy there. On the other hand, great diver- 
gence in the static temperatures is noted in the downstream region be y°”- 
the sphere-cone junction. This is due primarily to the shock shape which 
takes on an angle close to the respective semi- cone angles ( U c ). 1 he 

shock standoff distance normalized with respect to its own nose radius, 
i e A / R , also displays the effect of the semi- cone angle of the two 
bodies considered. However, T sh and A/R„ for the hemisphere -20° cone 
were essentially the same as a function of £ for both 280 and 310 kft. 


The temperature profiles within the viscous shock layer are shown 
in Figs. 36-38 for s 20* and. =9°. respectively. In the stagnation 

region, the maximum temperature is located immediately behind the shock, 
while in the downstream region the peak temperature moves toward the 
body surface, signifying the appearance of the viscous heating effect m a 
layer near the body. In a region sufficiently far downstream, this layer is 
thin compared to the shock-layer thickness, thus "recovering" a thin boundary- 
layer. Similar results were observed by Cheng (Ref. 15) and Davis (Ref. 70) in 
their analyses of hyperboloidal bodies. Thus a continuous solution is o tame 
which encompasses the "merged-layer " regime in the stagnation region through a 
"vo rticity -inte r action" regime (Refs. 16, 18) to a "boundary-layer regime in t e 
far downstream region of a blunted cone. Comparison of Figs. 36 and 37 with 

31 



Fig. 38 indicates a significant difference in the magnitude of T/T 

for a given value of g . This is primarily due to the difference in T Sh 

for the two bodies as illustrated on Fig. 35. 

The heat-transfer characteristics for a sphere cone at high altitudes 
are presented in Figure 39 in terms of the Stanton number (St ) along the 
body surface. The figure shows diminishing local heat-transfer rates in the 
downstream direction, as in Refs. 2 and 15. In addition, the results 
display the effect of the cone angle on the Stanton number such that the 
decrease is greater for the 0 C = 20° case than for the 0 C = 9° case. The 
influence of altitude can also be seen by comparing the O c - 20° results at 
310, 000 ft with those at 280, 000 ft. 

The distributions of the peak levels of the dissociated chemical species 
along the body are shown in Figs. 40 and 41 for the ^ = 20° cases. The rela- 
tive mass fractions of N£ and are not shown. Similar results were 
obtained for O c = 9° but are not included here. It is interesting to note from 
Figs. 40 and 41 that the present chemical nonequilibrium distributions for the 
neutral species (and also the ionized species) display a minimum far down- 
s tream (on the order of ten nose radii) of the junction point and show a gradual 
increase in their levels beyond this point. This type of behavior was observed 
for the case of a thin boundary- layer flow over a blunted cone at low altitudes, 
e. g. , Refs. 71, 72. The present results for a viscous shock layer appear to 
display similar bluntness effects on the flow. It is also important to note 
the relative magnitudes of the atomic ions and and the diatomic ion 

N-, + in the nose region. In addition, the O? ion is shown on Fig. 41 to be 
important in the downstream region for the 310 kft result's. These calculations 
suggest that all of the air ions should be included in the chemical model. 

The electron-density profiles calculated for a nonequilibrium, reacting 
gas are shown in Figs. 42 and 43 for Q c - 20° cases at various downstream 
locations. For the 280,000 ft calculation, the electron-density levels sharply 
decrease in the downstream direction until at about ten nose radii a 
minimum is reached due to expansion around the sphere- cone junction, and 
thereafter the number density gradually increases as mentioned previously. 

In the case of the 310, 000 ft calculations, this minimum occurred at about 
30 nose radii downstream. Similar results were also obtained for the 
& c ~ 9° case for both the nose region and the far downstream region but 
are not included here. 


5.4 Estimated Plane-Wave Transmission Coefficients for an 
Approximate Space-Shuttle Flow Field 

In order to obtain some estimate of the potential communications 
difficulties with the space- shuttle orbiter during earth entry, plane- wave 
transmission coefficients were calculated for several axial locations along 
the windward generator at several trajectory points. These calculations 
are very approximate at the present time because an accurate flow-field 
description for the orbiter is not available. However, by approximating 
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the flow field, the electron-density and collision-frequency distributions 
necessary for the E-M wave calculations were obtained, including non- 
equilibrium chemical reactions. Only pure air chemistry was included m 
these calculations and the chemical reactions and associated reaction rate 
coefficients are given in Table 1. The possible influence of ablation-product 
ionization on the electron-density level was not included. This consideration 
is particularly important at the lower altitudes where the electron density 
from pure air ionization is decreasing. 


In an effort to make maximum use of published work, the approximate 
flow field obtained by Lordi, Vidal and Johnson (Ref. 68) for the windward plane 
of symmetry of a blunted delta orbiter was used at altitudes of 250, 000 ft 
and below. These authors located streamlines in the inviscid flow field and 
determined the pressure distributions for these streamlines. They also 
performed calculations of the nonequilibrium species distributions at several 
altitudes of interest in this study for several of the streamlines. Because 
flow-field calculations were required over a large altitude range, we have 
used their calculations where available and performed additional nonequi- 
librium calculations at other trajectory points of interest by utilizing the 
streamline locations and pressure distributions supplied by Lordi and Vidal. 

It is important to note that the only ion included in their calculations was 
NO+, For their purposes, including NO + was appropriate but for the purposes 
of estimating electron-density levels, future calculations should include 
the remaining ions of air. Therefore, to be consistent, we used only the JNU 

ion in our low-altitude streamtube calculations. However, as previously 

. . r a_i tvt' O t M 1 anrl O . were 

U 2 » 


noted, the additional ions of the air model, N 2 , 
used in our high-altitude viscous-layer analysis. 


N , and O , were 


The plane- wave transmission coefficients as a function of altitude are 
shown in Fig. 44 to 48 for the windward generator at S/Rj^ = 1.22, 10, 30, 

50 and 90. Carrier frequencies of 100, 2200 and 10, 000 MHz have been con- 
sidered in this study. For altitudes of 250, 000 ft and below, the 20“ angle 
of attack flow-field description of Lordi, et. al. was used to r ® q ^ re 

electron-density and collision-frequency distributions. Above 250, 000 ** ttle 
viscous flow-field analysis for a 20° blunted cone described earlier in this 
report was used. For the purposes of the E-M wave claculations, the gradi- 
ent of electron density in the direction of propagation was accounted for by 
assuming the plasma layer to be composed of many smaller homogeneous 
slabs. However, the calculation technique (Ref. 73) is restricted in that it 
requires that the gradients normal to the direction of electromagnetic propagation 
be negligible . 


A typical shuttle entry trajectory and the trajectory points for which 
calculations were performed are given on Fig. 49- The calculations were 
based on an early trajectory and do not correspond exactly to the updated 
trajectory given on Fig. 49 • However, for the preliminary nature of this 
study, the correspondence is sufficiently close. It is significant to note that 
for an operating frequency of 2200 MHz, the transmission from an antenna 
located at X/R N of 1.22 will probably be significantly attenuated for a 
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period of approximately 15 minutes during the entry trajectory as a result 
of pure air ionization. 

The preliminary nature of the results described here must be empha- 
sized. They are presented here only for the purpose of estimating the severity 
of a potential problem. Before the selection of antenna frequencies and lo- 
cations is made, it is recommended that a more detailed description of the 
flow field should be obtained and that the potential influences of ablation product 
ionization should be assessed. 

As mentioned above, the results of the previously published flow field 
work of Lordi, Vidal and Johnson for the windward plane of symmetry of a 
blunted delta orbiter was used to obtain streamline locations and pressure 
distributions, as well as electron-density distributions for some streamlines, 
at altitudes below 250,000 ft. In order to provide a more complete report, 
the electron-density distributions obtained from these calculations are given 
in Figs. 50-53. The reader is cautioned that the only positive ion used in 
these calculations was NO+ and that ablation-product contamination at low 
altitudes was not considered. 


5. 5 Conclusions 

A theoretical analysis is presented which yields a continuous solution 
for a nonequilibrium viscous shock-layer flow around a sphere cone at high 
altitudes. Results are given for two different bodies ( © c = 9° and 20°) in 
terms of the temperature, the Stanton number and the electron-density distribu- 
tions. These calculations show the effects of the bluntness, chemical non- 
equilibrium and the semi- cone angle on the above quantities. The present 
analysis was undertaken as a first step in the possible application of an 
approximate technique to the calculation of a complicated flow field surrounding 
the space shuttle. 

Plane-wave transmission coefficients were calculated for several axial 
locations along the windward generator of a 20° half-angle spherically blunted 
cone and a blunted delta at 20° angle of attack. These calculations considered 
only pure air ionization. Transmission difficulties would probably be encounterei 
from approximately 300 to 200 kilofeet for frequencies less than 2200 MHz. 

Ab lation product ionization would tend to increase this altitude range. The 
transmission-coefficient results presented here are very approximate. Before 
the antenna frequencies and locations are selected for the shuttle vehicle, a 
more detailed analysis should be performed. 
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6 . POST -RETRACTION ELECTROSTATIC-PROBE RESULTS 


Voltaee-swept thin-wire electrostatic probes were flown on the RAM C-III 

vehicle as described in detail in Section 2. At an altitude of 'voltage 

200 000 ft the probes were retracted into the base-flow g & 

sweeping was continued. It was therefore the purpose of this study to tnvesUgate 
these post-retraction data and determine if useful electron-density and electron- 
temperature results could be obtained. 

The complex nature of the flow field in the base-flow region precluded 
use of the ion-current portion of the probe characteristic for determining 
the electron density. However, the electron-current portion of the characteristic 
could possibly be used. As discussed in Section 4, the collected current in 
this portion of the characteristic is independent of the ion mass and also inde- 
pendent of the probe orientation relative to the mass motionof heavyparUc le s . 
Because of the relative insensitivity of electron collection to the mag t 
of the iL-neutra! and neutral- neutral mean free paths, it is anticipated that 
the electron temperature deduced from the electron-retarding 
the electron current is much larger than the ion current, should be correct. 

The flight data were reviewed in some detail and it was determined 
that those data^n the altitude range of 190. 000 to 180, 000 ft we^ r 
free of foreign materials which were injected into the flow field as part of 
*rfli K ht-tes g t program. Therefore, it was felt that the probe charactertsttcs 
obtained in thi? alfitude range could be used to infer the 

and electron density in the base region. The technique used to obtain these 
values from the probe characteristics has been previously described in 
Sections 2 and 4 and will not be repeated here. 

Figures 54 and 55 illustrate the base region e \ ect / on - denS j7/^ o e ^ Ctr ° n - 
temperature data for the RAM C-III flight for the altitude range 190, 000 to 
180, P 000 ft. One of the four thin-wire probes was not operating during the 
Dost- retraction period so that results could only be obtained for three loca 
l “ns in the base P region instead of four. The location of the probe relative 
to the vehicle shoulder is illustrated on the figures The el ®f 0 rO - ,^ 3 ^ 
at 7 to 10 cm. of the shoulder was on the order of 3 to 4 x 10 e /cm b 
fell off rapidly to approximately 2 x 1 0 9 at 13 cm. It is difficult to compare 
these flight re y sults P with predicted values because of the complexities involved 
in determining the flow field in the near-wake region. However, one can 
make comparisons with electron densities calculated ior t e s ou f er , 
to expansion into the base region in order to obtain an indication of whether 
or not the results are reasonable. Huber has calculated the peak electron 
number density in the plasma layer prior to expansion to be approximately 
2 x 10 11 e"/cm 3 for the altitude range of interest. Qualitatively, the flight 
data nearest the shoulder appear to be compatable with the calculated v lu . 


* These unpublished results were generously supplied to us by Mr. Paul Huber 
of NASA, Langley Field. 
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Figure 55 suggests that the electron temperature in the base region 
was reasonably uniform between 8 and 13 cm. at a value of about 4000 to 
5000° K. The peak heavy-particle translational temperature in the plasma 
layer at the shoulder prior to expansion calculated by Huber was approximately 
4000° K. Once again, the flight data are qualitatively compatable with the 
calculations suggesting that the electron temperature in the base region may 
be nearly equal to the heavy-particle temperature just prior to the expansion. 
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7 THREE-DIMENSIONAL, VISCOUS SHOCK LAYER FOR SPACE 

SHUTTLE-TYPE BODY 

Section 5 of this report described preliminary estimates of the high 
attitude heat- transfer and Electron-density distributions ^ 
for the windward plane of symmetry by approximating the shuttle vehicle 
a blunted cone at zero angle of attack. These calculations were used to 
estimate the electromagnetic- wave transmission coef ^^ 6 

stations for a typical entry trajectory. The purpose of the work describe 
in this section wEis to formulate a method for obtaining solutions of the viscous 
shock u£r for a shuttle-type body at angle o£ attack ,c .that the 
results presented in Section 5 can be improved and so that calculations can 
be performed for locations other than the windward generators. 

The viscous hypersonic flow around a space shuttle vehicle at high 
altitudes at angle of attack is described by the Navier-Stokes equations (Ref. 74) 
and other continuum conservation equations (Ref. 75). The coordinate system 
chosen is "body oriented", similar to the conventional boundary-layer co- 
ordinate system (Ref. 18). 

Based on experience with the axisymmetric viscous-flow analyses 
(Refs. 2,20,41,69), the following assumptions are made: 

1) a thin shock layer 

2) two-layer model of Cheng (Ref. 15) 

3) Negligible changes in the flow properties such as the velocity 
and total enthalpy due to the chemical reaction in the flow field 

4) constant Prandtl and Schmidt numbers 

5) binary diffusion due to concentration gradient only 

6) negligible nonequilibrium effects due to electronic and vibrational 
relaxation 

7) ambipolar diffusion for the electrons and ions 

These assumptions have been made in part or wholly m the previous analyses ^ 
and the rationale for these assumptions may be found, for example, Refs. 

BASIC EQUATIONS 

With the above assumptions the governing equations for the three- 
dimensional, viscous, ionized shock layer become s (see Figure 1 for a descrip- 
tion of the flow field). 
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Streamwise Momentum: 
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With the two- layer model of Cheng, these equations are applied to the 
viscous shock layer allowing for the diffusion of the species into the shock- 
transition zone by using the modified Rankine-Hugoniot conditions. The 
boundary conditions for a fully catalytic solid wall are: 
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T ransf ormation 

In solving the equations the following coordinate transformations are 
made: . « _ pL* . ^ n s r oi 4> 
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In addition two stream functions are introduced such that 
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Along with the nondimensional velocities 1 a/ defined as 
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With introduction of the above transformation, the equations become, after 
a series of manipulations: 
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Transverse Momentum: 
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At the conclusion of the present contract, the formulation has been 
completed but results have not been obtained. However, possible methods of 
solution have been investigated and some tentative conclusion regarding the 
most promising approaches are presented. 


As methods of solution to the partial, non-linear differential equations, 
the integral method (Refs. 2,20,21,69), the unsteady method (Refs . 76, 77) , the 
explicit method (Refs. 78,79), and the implicit method (Refs. 70,80) were con- 
sidered. From the viewpoint of quick but approximate results, the integral 
method has been shown (Refs. 2,20,21,69) to be adequate, even in the three- 
dimensional flow situation (Ref. 81). If, however, more accurate and 
elaborate results are desired, the implicit scheme appears to be the 
most reasonable method, being usually free from convergence problem 
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even in the non-linear case. For the case of a three-dimensional imcompres- 
sible, boundary- layer flow, initial success using the implicit scheme is 
reported in Ref. 80 and 82. However, because the geometry of the space 
shuttle is rather complicated, both the integral method and the implicit scheme 
should be investigated before set tling on one particular method. 
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APPENDIX A 


Admittance data taken during the experiments are presented relative 
to a short circuit as are the values derived from the circular aperture 
theoretical model. In special cases, where the antenna is experimentally 
matched to free space or partially matched using tuning stubs, calculations 
must be adjusted accordingly*. 

No adjustments are necessary for RAM C-III antenna since the 
circular aperture theoretical model values include the natural mismatch of 
the aperture and plug. However, because the RAM C-II antenna admittance 
for the no plasma case was different from its natural value, adjustments 
were necessary. The equivalent circuit used to make the adjustments is 
shown below. 


Turing Tuning- ft perturb 

Stub Stub Pubne 

#2 #1 



Figure A-l 


The aperture admittance Xx. is transformed to a line length p JL , 
to the first stub and the shunt susceptance, -^ s , added to the network. Both 

P£ 1 and have values such that the aperture admittance of the antenna 

is matched to free space in the no plasma case. The admittance of the antenna 
at the terminals of the first stub is given by the transmission line equation 


y*- + i tayn ' ft jj. + • £ 
i + i y«, * s 


* The authors are grateful to Dr. John Q. Howell for formulating these 
relationships and applying them to the unmatched antenna. 
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A second stub is added at a distance pfi^ ^ rom the fi rs ^ stub such 
that the transformed admittance at the measurement reference plane is the 
same value as measured in the experiment for the no plasma case. 

The position of the reference plane is the null location for a shorted 
aperture and is adjusted by selection of p-fc^ a fter the length PA 3 is deter- . 
mined, p £3 and 4* Sz were adjusted to give an admittance value. of 1.38 .18 

at the reference plane for the no plasma case. Only minor deviations from this 
value were observed during the course of the experiments as shown in Table II. 
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Table! 


CHEMICAL REACTIONS AND RATE COEFFICIENTS USED 
IN NONEQUILIBRIUM CALCULATIONS 


NO. 

REACTION 

FORWARD RATE COEFF, k F 

BACKWARD RATE COEFF, k B 

THIRD BODY, M 


FORWARD DIRECTION 

cm ^/m ole sec 

O e o 

cm /mole sec OR cm/mole sec 


1 

0 2 + M -► 20 + M 

3.6 x 10 18 T -10 expl— 5.95 x 10 4 /T) 

3.0 x 10 15 T -0 - 5 

N, NO 

2 

N 2 + M — 2N + M 

1.9 x 10 17 T -0,5 expl— 1.13 x 10 5 /T) 

1.1 x 10 16 T -0 ' 5 

O, NO, 0 2 

3 

NO + M N + 0 + M 

3.9 x 10 20 T -1,5 exp(-7.55 x 10 4 /T) 

1.0 x 10 20 T -1 ' 8 

z 

CS) 

O 

4 

O + NO — ► N + 0 2 

3.2 x 10 9 T 1 expl— 1.97 x 10 4 /T> 

1.3 x 10 10 T 10 expl— 3.58 x 10 3 /T) 


5 

0 + N 2 N + NO 

7.0 x 10 13 expl— 3.8 x IoVt) 

1.56 x 10 13 


6 

N + N 2 — ► N + N + N 

4.085 x 10 22 T -1 - 5 exp(— 1.13 x 10 5 /T) 

2.27 x 10 21 T -1 5 


7 

0 + N -► NO + + e“ 

(1.4 + 0.4) x 10 6 T 1 ' 5 exp(— 3.19 x 10 4 /T) 

(6.7 + 2.3) x 10 21 T -1 - 5 


8 

0 + e~ ♦ 0 + + e + e 

(3.6+ 1.2) x 10 31 T -2 * 91 exp(-1 58 x10 5 /T) 

(2.2 + 0.7) x 10 40 T^ 4,5 


9 

N + e"~ N + + e + e 

(1.1 + 0.4) x 10 32 T -3 - 14 expl— 1.69 x 10 5 /T) 

(2.2 + 0.7) x 10 40 T -4 ' 5 


10 

O + O -► 0 2 + + e” 

(1.6 + 0.4) x 10 17 T ' 0 98 expl— 8.08 x 10 4 /T) 

(8.0 + 2.0) x ID 21 T -1 - 5 


11 

o + o 2 + -► o 2 + 0 + 

2.92 x 10 18 T -1 ’ 11 expl— 2.8 x 10 4 /T) 

7.8 x 10 11 T 0 - 5 


12 

N 2 + N + -► N + N 2 + 

2.02 x 10 11 T 0 - 81 expl— 1.3 x 10 4 /T) 

7.8 x 10 11 T 0 - 5 


13 

N + N “► N 2 + + e” 

(1.4 + 0.3) x 10 13 exp(-6.78 x 10+/T) 

11.5 + 0.5) x 10 22 T -15 


14 

° 2 + N 2 NO + NO + + e — 

1.38 x 10 20 T -1 - 84 expl— 1.41 x 10 5 /D 

1.0 x 10 24 T -2 ' 5 


15 

NO + N 2 — ► NO + + e _ + N 2 

2.2 x 10 15 T -0 - 35 expl— 1.08 x 10 5 /T) 

2.2 x 10 26 T -2 5 


16 

0 + NO + — ► NO + 0 + 

3.63 x 10 15 T -0 - 6 expl— 5.08 x 10 4 /T) 

1.5 x 10 13 


17 

n 2 + 0 + -► o + n 2 + 

3.4 x 10 19 T -2 0 expl— 2.3 x 10 4 /T) 

2.48 x 10 19 T -2-2 


18 

N + NO + — ► NO + N + 

1.0 x 10 19 T -0 - 93 expl— 6.1 x 10 4 /T) 

4.8 x 10 14 


19 

0 2 + NO*-fr- NO + 0 2 + 

1.8 x 10 15 T 0 ’ 17 expl— 3.3 x 10 4 /T) 

1.8 x 10 13 T 0 - 5 


20 

0 + NO + -► 0 2 + N + 

1.34 x 10 13 T 0 - 31 expl -7.727 x 10 4 /T) 

1.0 x 10 14 


21 

NO + 0 2 — ► NO + + e~ + 0 2 

8.8 x 10 15 T -0 ' 35 expl -1.08 x 10 5 /T) 

8.8 x 10 26 T -2,5 


22 

0 2 + 0 — ► 20 + 0 

9.0 x 10 19 T -10 expl -5.95 x 10 4 /T) 

7.5 x 10 16 T -0 - 5 


23 

°2 + °2 20 + °2 

3.24 x 10 19 T -10 expl -5.95 x 10 4 /T) 

2.7 x 10 16 T -0 - 5 . 


24 

0 2 + N 2 -► 20 + N 2 

7.2 x 10 18 T — 10 expl— 5.95 x 10 4 /T) 

6.0 x 10 15 T -0 ' 5 


25 

N 2 + N 2 -► 2N + N 2 

4.7 x 10 17 T -0 - 5 exp(-1.13 x 10 S /T) 

2.72 x 10 16 T -0 - 5 


26 

NO + M — ► N + 0 + M 

7.8 x 10 20 T -1 - 5 expl— 7.55 x 10 4 /T) 

2.0 x 10 20 T -1 - 5 

O, N, NO 



Table II 

SUMMARY OF EXPERIMENTAL RESULTS OBTAINED WITH RAM C-IIIS-BAND ANTENNA 
AT 6.75 INCHES FROM PLATE LEADING EDGE IN CARBON MONOXIDE PLASMA 


RUN 

# 

PEAK n e OVER 
ANTENNA 
e'/cm 3 

n e /r VcR 

ANTENNA IN FREE SPACE 
BEFORE PLASMA 

ANTENNA IN FREE 
SPACE ROTATED & 
CORRECTED FOR 
LINE LOSS 

ANTENNA WITH P 
OVER APERTU 

LASMA 

RE 


ANTENNA WITH 
PLASMA ROTATED & 
nnRRECTED FOR 

#1 

ARC 

#2 

ARC 

#3 

ARC 

#4 

ARC 

REFL. 

ARC 

ARC 

#2 

ARC 

#3 

ARC 

#4 

ARC 

REFL. 

ARC 

LINE LOS 

S 

Y 

in 

Y 

in 

FOIL 

UNDER PLEXIGLAS 

L — 

1.872 

1.08 

0.342 

1.58 

0.90 










15 

(6.3 TO 5.3) x 10 11 

4.2 + .3 

0.716 

1.14 

1.30 

0.944 

0.318 

0.48 + j 0.02 

0.35 

0.476 

1.472 

1.764 

1.070 

0.76 

0.072 + j 0.05 

0.862 

16 

(7.6 TO 6.4) x 10 11 

5.0 + .4 

0.728 

1.164 

1.32 

0.952 

0.325 

0.47 - 0.001 

0.36 

0.610 

1.570 

1.74 

0.956 

0.772 

0.077 - j 0.09 

0.860 











0.464 

1.476 

1.716 

0.896 

0.762 

0.092 - j 0.038 

0.830 


*(5.4 TO 4.8) x 10’ 

3.6 + .3 















17 

*(4.0 TO 3.4) x 10 11 

2.6 + .2 

0.734 

1.144 

1.306 

0.956 

0.31 

0.49 + j 0.005 

0.342 

0.376 

1.404 

1.816 

1.230 

0.80 

0.060 + j 0.085 

0.890 


(3.6 TO 3.0) x 10 11 

2.4 + .2 i 








0.324 

1.304 

1.780 

1.290 

0.792 

0.065+ j 0.147 

0.880 


‘(4.0 TO 3.4) x 10 11 

2.6 + .2 








0.276 

1.396 

1.764 

1.142 

0.792 

0.063+ j 0.06 

0.880 

18 

(8.2 TO 7.0) x 10 11 

5.4 + .45 

0.722 

1.12 

1.284 

0.944 

0.302 

0.50 + j 0.025 

0.33 

0.590 

1.520 

1.720 

0.974 

0.760 

0.085 - j 0.075 

0.845 

19 

(1.3 TO 1.2) x 10 11 

0.9 + .15 

0.704 

1.136 

1.314 

0.958 

0.320 

0.48 + j 0.03 

0.35 

0.514 

1.040 

1.674 

1.444 

0.710 

0.14 +j 0.418 

0.790 


*(1.05 TO 1.0) x 10 11 

0.74 + .12 








0.682 

0.878 

1.444 

1.380 

0.525 

0.34 + j 0 .49 

0.582 

20 

(9.0 TO 7.5) x 10 10 

0.6 + .1 

0.704 

1.136 

1.310 

0.96 

0.318 

0.48 + j 0.03 

0.35 

0.796 

1.020 

1.240 

1.056 

0.249 

0.625 + j 0.15 

0.25 


NOTE: 1) THE VOLTAGE REFLECTION COEFFICIENT, \P\, GIVEN HERE HAS BEEN CORRECTED FOR LINE LOSS. TO GET THE UNCORRECTED VALUE, MULTIPLY BY 0.90. 

2 ) * INDICATES POST TEST-TIME DATA AND t INDICATES PRE TEST-TIME DATA. 




CJI 

ft- 


Table HI 

SUMMARY OF EXPERIMENTAL RESULTS OBTAINED WITH RAM C-IL S-BAND ANTENNA 
AT 11.25 INCHES FROM PLATE LEADING EDGE IN CARBON MONOXIDE PLASMA 


RUN 

# 

PEAK 

ELECTRON 

DENSITY 

OVER 

ANTENNA 

e'/cm’ 3 

n e /n e^CR 

ANTENNA IN FREE SPACE BEFORE 
PLASMA ARRIVAL 

ANTENNA IN FREE SPACE 
ROTATED AND CORRECTED 
FOR LINE LOSS 

ANTENNA WITH PLASMA 
OVER APERTURE 


ANTENNA WITH PLASMA 
ROTATED AND 
CORRECTED FOR LINE LOSS 

#1 

ARC 

#2 

ARC 

#3 

ARC 

#4 

ARC 

REFL. 

ARC 

#1 

ARC 

#2 

ARC 

#3 

ARC 

#4 

ARC 

REEL. 

Y 

in 

ARC 

Y 

in 







SHORT 

OVER ANTENNA AF 

ERTURE 

1.280 

0.180 

1.325 

1.830 

0.843 










1 

(2.46 + .24) x 10 11 

1.76 + .17 

0.95 

0.89 

1.07 

1.12 

0.148 

1.35 + j 0.19 

0.17 

1.376 

1.498 

0.916 

0.630 

0.565 

0.205 - j 0.246 

0.68 

2 

(2.87 + .28) x 10 11 

2.05 + .20 

0.96 

0.864 

1.078 

1.14 

0.148 

1.38 + j 0.18 

0.18 

1.522 

1.396 

0.670 

0.842 

0.630 

0.22 - j 0.470 

0.705 

3 

(1.46 + .13) x 10 11 

1.04 + .09 

0.964 

0.874 

1.078 

1.152 

0.140 

1.39 + j 0.19 

0.18 

0.722 

1.516 

1.572 

0.774 

0.600 

0.17 + j 0.450 

0.760 

4 

<1.60+ .12) x 10 11 

1.14 + .09 

0.976 

0.890 

1.086 

0.144 

0.131 

1.36 + j 0.18 

0.165 

1.136 

1.620 

1.200 

0.434 

0.58 

0.175 + j 0.04 

0.710 

* 

(1.13 + .09) x 10 11 

0.81 + .06 

0.976 

0.890 

1.086 

0.144 

0.131 

1.36 + j 0.18 

0.165 

0.58 

1.124 

1.550 

1.176 

0.515 

0.49+ j 0.97 

0.617 

5 

(1.60 + .12) x 10 11 

1.14 + .09 

0.970 

0.862 

1.050 

1.134 

0.130 

1.38 + ] 0.12 

0.165 

1.196 

1.550 

1.170 

0.476 

0.56 

0.21 - j 0.016 

0.652 

* 

(1.03 + .05) x 10 11 

0.73 + .03 

0.970 

0.862 

1.050 

1.134 1 

0.130 

1.38 + j 0.12 

0.165 

0.88 

0.854 

1.240 

1.230 

0.360 

1.72 + j 0.85 

0.39 

6 

(4.07 + .73) x 10 10 

0.29 + .05 

0.976 

0.876 

1.050 

1.140 

0.122 

1.38 + j 0.10 

0.165 

0.976 

0.876 

1.050 

1.140 

0.122 

1.38 + j 0.10 

0.165 

7 

(3.11 + .3) x 10 10 

0.22 + .02 

0.944 

0.876 

1.064 

1.122 

0.140 

1.33 + j 0.24 

0.18 

0.996 

0.904 

1.090 

1.120 

0.140 

1.30+ j 0.13 

0.145 

8 

(4.29 + .51) x 10 10 

0.31 + .03 

0.944 

0.861 

1.058 

1.122 

0.145 

1.38 + j 0.18 

0.175 

0.944 

0.861 

1.058 

1.122 

0.145 

1.38 + j 0.18 

0.175 

9 

(7.48 + .82) x 10 10 

0.53 + .06 

0.944 

0.864 

1.050 

1.122 

0.143 

1.36 + j 0.17 

0.17 

0.830 

0.900 

1.232 

1.160 

0.280 

1.19 + j 0.70 

0.315 

* 

(6.62 + .72) x 10 10 

0.47 + .05 

0.944 

0.864 

1.050 

1.122 

0.143 

1.36 + j 0.17 

0.17 

0.856 

0.762 

1.110 

1.240 

0.305 

1.70 + j 0.65 

0.341 

10 

(7.68 + .73) x 10 10 

0.55 + .05 

0.976 

0.876 

1.050 

1.136 

0.142 

1.38 x j 0.10 

0.165 

0.856 

0.840 

1.312 

1.296 

0.415 

1.51 + j0.95 

0.40 

11 

(1.30 + .12) x 10 11 

0.93 + .08 

0.966 

0.836 

1.038 

1.136 

. 0.141 

1.43 + i 0.10 

0.182 

0.57 

1.336 

1.562 

0.952 

0.560 

0.26 + i 0.67 

0.70 


NOTE: 1) THE VOLTAGE REFLECTION COEFFICIENT, (PI, GIVEN HERE HAS BEEN CORRECTED FOR LINE LOSS TO GET THE UNCORRECTED VALUE, MULTIPLY BY 0.835. 

2) ANTENNA FREQUENCY WAS 3348 MHz AND THEREFORE n e l CR ^1.4 10 1 1 e'/cm 3 

3) * INDICATES POST TEST-TIME DATA 




Table GZ 

RESULTS OBTAINED FROM LANGMUIR PROBE CURRENT-VOLTAGE CHARACTERISTICS 


TEST 


M -WAVE NUMBER 

NUMBER DENSITY FROM 


T e 

NUMBER DENSITY FROM 

V X D 

< v oS v fW 

(V ooV exp . 

GAS 

PROBE LOCATION 

DENSITY AT 22.5 in. 

ION SATURATION 

Rp/Xn 


ELECTRON SATURATION 

- 

- 

- 

e~ /cm** 

e /cm^ 

ol 

y 

e'/cm^ 

— 

volts 

volts 

CO 

FREE STREAM 

2.0 x 10 11 

1.5 x 10 11 

7.5 

2110 

1.3 x 10 11 

6.9 

0.9 

■g 

CO 



1.5 x 10 11 

1.7 x 10 11 

8.8 

1650 

1.7 x 10 11 

8.7 

0.7 


CO 

T 

BOUNDA 


2.0 x 10 11 

1.2 x 10 11 

7.0 

1945 

1.4 x 10 11 

7.4 

0.85 

0.8 

CO 

RY LAYER 


1.4 x 10 10 

1.7 


1.3 x 10 10 

1.7 

1.2 

1.2 

CO 




1.3 x 10 10 

1.9 


1.3 x 10 10 

1.8 

0.9 

1.1 

CO 




1.7 x 10 10 

2.0 

Sv '-Til 

1.4 x 10 10 

1.4 

1.0 

1.1 

CO 




4.4 x 10 10 

3.3 


3.3 x 10 10 

2.8 

1.0 

1.2 

CO 




1.6 x 10 10 

1.9 

2510 

1.4 x 10 10 

1.8 

1.1 * 

0.9 

CO 




6.6 x 1C 10 

4.3 

E31 

5.3 x 10 10 

3.8 


0.9 

CO 




1.1 x 10 10 

1.6 

2170 

1.1 x 10 10 

1.7 

0.95 

0.9 

CO 




5.3 x 10 10 

3.5 

2310 

5.0 x 10 10 

3.4 

1.0 

1.2 

CO 




1.8 x 10 10 

2.3 

1810 

1.8 x 10 10 

2.3 

0.8 

0.8 

CO 




6.2 x 10 10 

4.4 

1730 

6.5 x 10 10 

3.6 

0.76 

0.8 

CO 




Z1 x 10 10 

2.3 

2100 

2.1 x 10 10 

2.3 

0.92 

1.0 

CO 




4.1 x 10 10 

3.4 

1950 

4.1 x 10 10 

3.4 

sa 

0.9 

CO 




8.1 x 10 9 

1.5 

1955 

8.9 x 10® 

1.6 

1 

0.8 

CO 




5.0 x 10 10 

3.9 

1960 

4.6 x 10 10 

3.6 

ES 


CO 




2.5 x 10 10 

2.5 

2130 

2.6 x 10 10 

2.5 

El 

1.1 

CO 




9.4 x 10 9 

1.6 


9.2 x 10® 

1.6 

Bni 

0.8 

CO 




2.7 x 10 10 

2.7 

2060 

2.9 x 10 10 

2.8 

19 

1.0 

CO 




1.2 x 10 10 

1.7 

2160 

1.3 x 10 10 

1.8 

$t:3 

0.9 

CO 




2.8 x 10 10 

2.8 

1980 

3.8 x 10 10 

3.2 

0.87 

1.0 

CO 




6.8 x 10® 

1.3 

2160 

6.5 x 10® 

1.3 

0.95 

1.1 

CO 


r 


4.3 x 10 10 

3.5 

1915 

4.3 x 10 10 

1.4 

0.84 

1.0 

N 2 

F 

FREE STREAM 

9.5 x 10 9 

9.0 x 10® 

1.4 

3580 

1.1 x 10 10 

1.5 

1.6 

2.0 

N 2 



1.0 x 10 10 

1.1 x 10 10 

1.5 

3670 

1.1 x 10 10 

1.5 

1.6 

1.6 

N 2 



15 x 10 10 

14 x 10 10 

1.7 

3510 

1.1 x 10 10 

1.5 

1.53 

1.5 

n 2 


r 

7.2 x 10 9 

7.4 x 10® 

1.5 

2450 

3.5 x 10® 

2.1 

1.06 

0.7 

















Figure 1 SCHEMATIC FLOW FIELD 



Figure 2 RAM C-E PAYLOAD EXTERNAL 
CONFIGURATION 
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1 

0.378 

(0.96) 

2 

0.716 

(1.84) 

3 

1.062 

(2.70) 

4 

1.408 

(3.58) 

5 

1.754 

(4.46) 

6 

2.102 

(5.34) 

7 

2.438 

(6.20) 

8 

2.790 

(7.08) 


Figure 3 RAM C-H CONSTANT BIAS- 
VOLTAGE ELECTROSTATIC 
PROBE RAKE 
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PHENOLIC —GRAPHITE 
NOSE COVERING 
6,28 in. (IS, 05 cm) radius 


TEFLON HEAT; 
SHIELD ON 
AFTERBODY 


i CONSTANT BIAS'- 
VOLTAGE PROBE 
RAKE LOCATED 
AT X/O - 3.87 


SWEPT BIAS- 
VOLTAGE PROBE 
RAKE LOCATED 
| AT X/P 3.87 


Figure 4 RAM C-IH PAYLOAD EXTERNAL 
CONFIGURATION 
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BERYLLIUM OXIDE 
INSULATOR 


-0.005 (0.0127) radius 


PHENOLIC 

FIBERGLASS 

HOLDER 

k 



PAYLOAD 

SURFACE 


ION 

COLLECTOR 

STANDOFF DISTANCE, y 
(in.) (cm) 

1 

0.35 

0.88 

2 

0.69 

1.75 

3 

1.04 


4 

1.38 


5 

1.72 


6 

2.07 


7 

2.42 


8 

2.76 


9 

3.10 


10 

3.45 


11 

3.80 


12 

4.14 


13 

4.48 


14 

4.83 


15 

5.18 

13.15 

16 

5.52 

14.02 


Figure 5 RAM C-IH CONSTANT BIAS- 
VOLTAGE ELECTROSTATIC 
PROBE RAKE 
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-0.005 (0.0127) diameter 
0.400 (1.016) length 

0.040 (0.1016) 



BERYLLIUM 

OXIDE 

INSULATOR 


PHENOLIC 

FIBERGLASS 

HOLDER 


LEADING-EDGE 

DETAIL 


PAYLOAD 

SURFACE 



ELECTRON 

COLLECTOR 

STANDOFF DISTANCE, y 
(in.) (cm) 

1 

0.59 

1.5 

2 

1.38 

3.5 

3 

2.56 

6.5 

4 

3.74 

9.5 


Figure 6 RAM C-HI SWEPT-VOLTAGE 
ELECTROSTATIC PROBE RAKE 
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NON-EQUILIBRIUM FLOW, 
VISCOUS-SHOCK LAYER THEORY 

KEY TO FLIGHT DATA 


VsNS^ RAM C-HI CONSTANT BIAS- 
VOLTAGE PROBE 


T RAM C-IC CONSTANT BIAS- 
1 VOLTAGE PROBE 



Figure 7 COMPARISON BETWEEN CAL- 
CULATED ELECTRON DENSITY 
AND IN-FLIGHT ELECTROSTATIC 
PROBE MEASUREMENTS 
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TEMPERATURE 



280 270 260 250 240 230 220 210 

ALTITUDE, kilofeet 

Figure 8 COMPARISON BETWEEN CALCULATED 
HEAVY-PARTICLE TRANSLATIONAL 
TEMPERATURE AND IN-FLIGHT 
MEASURED ELECTRON TEMPERATURE 
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KEY TO VISCOUS SHOCK-LAYER THEORY 


UPPER -BOUND REACTION RATE COEFFICIENTS 
FOR DE-IONIZATION REACTIONS 


LOWER-BOUND REACTION RATE COEFFICIENTS 

FOR DE-IONIZATION REACTIONS 

KEY TO FLIGHT DATA 



Figure 9 INFLUENCE OF REACTION RATE 
COEFFICIENTS ON CALCULATED 
ELECTRON DENSITY AT ELECTRO- 
STATIC PROBE LOCATION 
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ELECTRON NUMBER DENSITY, e"/cm 


(SHOCK) 



Figure 10 INFLUENCE OF CHEMICAL MODEL 
ON CALCULATED ELECTRON 
DENSITY AT ELECTROSTATIC 
PROBE LOCATION 
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NOTE: RESULTS ARE FOR UPPER-BOUND 
REACTION RATE COEFFICIENTS 
FOR DE-IONIZATION REACTIONS 



Figure 11 NONEQUILIBRIUM SPECIES DIS- 
TRIBUTIONS IN PLASMA LAYER 
FOR 233,000 ft ALTITUDE AT 
25,000 ft/sec VELOCITY 
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PROBE LOCATION 



NOTE: RESULTS ARE FOR UPPER-BOUND 
REACTION RATE COEFFICIENTS 
FOR DE-IONIZATION REACTIONS 



Figure 12 NONEQUILIBRIUM SPECIES DIS- 
TRIBUTIONS IN PLASMA LAYER 
FOR 275,000 ft ALTITUDE AT 
25,000 ft/sec VELOCITY 
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ELECTRON NUMBER DENSITY, 


NON-EQUILIBRIUM FLOW, 
VISCOUS SHOCK-LAYER THEORY 


KEY TO FLIGHT DATA 
X'SS RAM C-m CONSTANT BIAS-VOLTAGE 
PROBE 

T RAM C-It CONSTANT BIAS-VOLTAGE 
1 PROBE 



DISTANCE NORMAL TO BODY, cm 

Figure 13 COMPARISON BETWEEN CAL- 
CULATED ELECTRON DENSITY 
AND IN-FLIGHT ELECTROSTATIC 
PROBE MEASUREMENTS AT 
214,000 ft ALTITUDE 
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ELECTRON NUMBER DENSITY, e/cm 



BODY DISTANCE NORMAL TO BODY/SHOCK- 
LAYER THICKNESS, y/A 

Figure 14 COMPARISON OF PREDICTED 
ELECTRON DENSITY PROFILE 
WITH PEAK VALUE FROM 
REFLECTOMETER FLIGHT 
DATA AT X/R n = 5.2 
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LAYER THICKNESS, y/A 


Figure 15 COMPARISON OF PREDICTED 
ELECTRON-DENSITY PROFILE 
WITH PEAK VALUE FROM 
REFLECTOMETER FLIGHT 
DATA AT X/R n = 7.65 
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PLEXIGLAS LAYER 



Figure 16 SKETCH OF THE RAM C-UI S-BAND 
ANTENNA 
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Figure 17 ASSUMED ELECTRON-DENSITY 

PROFILES FOR PRE-EXPERIMENT 
CALCULATIONS 
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Figure 18 INFLUENCE OF BOUNDARY LAYER 
ELECTRON - DENSITY PROFILE ON 
CALCULATED ADMITTANCE 
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Figure 19 (a)RAM CHI AND RAM C-H S-BAND 

AND X-BAND ANTENNAS MOUNTED 
IN FLAT PLATE 



Figure 19 (b) RAM C-IH AND RAM C-H ANTENNAS 
IN FLAT-PLATE MODEL 
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INCIDENT 




Figure 20 SCHEMATIC DIAGRAM OF FOUR- 
PROBE REFLECTOMETER 
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Figure 21 CALCULATED MEAN FREE PATHS IN 
CARBON MONOXIDE PLASMA 
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Figure 22 ELECTRON NUMBER DENSITY ABOVE 
PLATE AT 6.75-INCHES FROM 
LEADING EDGE IN CARBON 
MONOXIDE PLASMA 
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10 3 10 4 
ELECTRON TEMPERATURE, °K 


Figure 23 ELECTRON TEMPERATURE ABOVE 
PLATE AT 6.75-INCHES FROM 
LEADING EDGE IN CARBON 
MONOXIDE PLASMA 
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Figure 24 ELECTRON NUMBER DENSITY ABOVE 
PLATE AT 11.25 INCHES FROM 
LEADING EDGE IN CARBON 
MONOXIDE PLASMA 
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Figure 25 ELECTRON-TEMPERATURE ABOVE 
PLATE AT 11.25-INCHES FROM 
LEADING EDGE IN CARBON 
MONOXIDE PLASMA 
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(a) REFLECTOMETER DIODE #1 < b > REFLECTOMETER DIODE #2 
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(c) REFLECTOMETER DIODE #3 



(d) REFLECTOMETER DIODE #4 


1) ANTENNA LOCATED AT 6.75 in. FROM PLATE LEADING EDGE 

2) PEAK ELECTRON DENSITY OVER ANTENNA DURING TEST TIME WAS n e % 7.6 x 10 e*/cm 

3) RUN NO. 16, TABLE II 


Figure 26 TYPICAL REFLECTOMETER AND MICROWAVE-INTERFEROMETER OSCILLOSCOPE 
RECORDS FOR RAM C-III S-BAND ANTENNA IN CARBON MONOXIDE PLASMA 








(e) REFLECTOMETER REFLECTED SIGNAL 


Figure 26 TYPICAL REFLECTOMETER AND 
MICROWAVE-INTERFEROMETER 
OSCILLOSCOPE RECORDS FOR 
RAM C-HLS-BAND ANTENNA IN 
CARBON MONOXIDE PLASMA (Cont.) 
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I SMITH CHART FORM 42B3PH (2-4») 1 RAY ELECTRIC 


COMPANY. PINE BROOK. N.J. OlM» P AINTEP IN ULSA. | 
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I H~ THEORY FOR RAM C-HE ANTENNA 
FOR ■?/& = 0.1 AND n 

"peak 

AT y = 1.2 inche* FROM SURFACE 

('*0*') EXPERIMENTAL DATA .// 

OBTAINED FROM 4-PROBE (/ 

REFLECTOMETER x 
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Figure 27 COMPARISON OF EXPERIMENTAL 

DATA WITH THEORY FOR RAM C-HL 
ANTENNA AT 6.75 INCHES FROM 
PLATE LEADING EDGE 




Figure 28 COMPARISON OF EXPERIMENTAL 

DATA WITH THEORY FOR RAM C-H 
ANTENNA AT 11.25 INCHES FROM 
PLATE LEADING EDGE 
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DIMENSIONLESS CURRENT, 



(Rp/ X q) 2 i+ 

Figure 29 LAFRAM BOISE THEORETICAL RESULT 
FOR ION CURRENT COLLECTION 
WHEN Tj/T e = 0 
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(Rp/ X d) 2 *- 


Figure 30 LAFRAMBOISE THEORETICAL RESULT 
FOR ELECTRON COLLECTION 
WHEN Tj/T e = 0 
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PROBE CURRENT, ma 


O EXPERIMENTAL DATA 

LINE FAIRED THROUGH DATA 

RESERVOIR TEMPERATURE: 7060°K 
LOCATION: IN BOUNDARY LAYER, 11.25 in. FROM 
PLATE LEADING EDGE 



Figure 31 ION-CURRENT REGION OF 
EXPERIMENTAL PROBE 
CHARACTERISTIC IN 
CARBON MONOXIDE 
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PROBE CURRENT, ma 
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EXPERIMENTAL DATA 
LINE FAIRED THROUGH DATA 

RESERVOIR TEMPERATURE: 7060°K 
RESERVOIR PRESSURE: 17.3 atm 
LOCATION: IN BOUNDARY LAYER. 

11.25 in. FROM PLATE 
LEADING EDGE 
FROM EXPERIMENTAL DATA: 


= 5.37 
= 2160°K 

= 0.50 volts 


- 0.4 
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2.0 2.4 


PROBE VOLTAGE, volts 


Figure 32 ELECTRON-RETARDING AND 

ELECTRON-CURRENT REGIONS 
OF EXPERIMENTAL PROBE 
CHARACTERISTICS IN CARBON 
MONOXIDE 
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NORMALIZED CURRENT, i 



DIMENSIONLESS POTENTIAL, 
X p = (e/k T e ) (Vp - ) 


Figure 33 COMPARISON OF ION-CURRENT AND 
ELECTRON-RETARDING REGIONS 
WITH LAFRAMBOISE'S THEORY 
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DIMENSIONLESS POTENTIAL, 
Xp = (e/k T e ) (V p - ) 


Figure 34 COMPARISON OF ELECTRON- 
RETARDING AND ELECTRON- 
CURRENT REGIONS WITH 
LAFRAMBOISE'S THEORY 
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RE BEHIND SHOCK, T 



0 1 2 10 20 30 40 50 60 70 80 90 

NONDIMENSIONAL DISTANCE ALONG BODY, 


£ = x/R n 

Figure 35 TEMPERATURE BEHIND SHOCK 
AND SHOCK-LAYER THICKNESS 
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ALTITUDE = 280,000 ft 
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ALTITUDE = 233,000 ft 

= 25,100 ft/sec 



BEHIND SHOCK, T/Tg H 

Figure 38 TEMPERATURE PROFILE FOR 
HEMISPHERE - 9° CONE 
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NONDIMENSIONAL DISTANCE ALONG BODY, 
£ = x/R n 


Figure 39 STANTON-NUMBER DISTRIBUTION 
FOR HEMISPHERE-CONE BODY 
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MAXIMUM SPECIES MASS FRACTION, C: 

gmj/gm mix max 


(0 2 , N 2 NOT SHOWN) 



J =x/R n 


Figure 40 NONEQUILIBRIUM FLOW SPECIES 

DISTRIBUTION FOR HEMISPHERE - 
20° CONE 
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MAXIMUM SPECIES MASS FRACTION, C { . gm-Jgm mix 



NONDIMENSIONAL DISTANCE ALONG BODY, £ = x/R N 


Figure 41 NONEQUILIBRIUM FLOW SPECIES DISTRIBUTION FOR HEMISPHERE - 20° CONE 
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ELECTRON NUMBER DENSITY, e7cm 



DISTANCE NORMAL TO BODY/SHOCK-LAYER 
THICKNESS, y/A 


Figure 42 ELECTRON NUMBER DENSITY 

DISTRIBUTION IN SHOCK-LAYER 
FOR HEMISPHERE - 20° CONE 
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$ - 1.2 (JUNCTION) 

ALTITUDE = 310,000 ft 

U^ = 25,100 ft/sec 
R n = 4.0 ft 

£ - 3 = x/R n 

SHOCK 


.2 0.3 0.4 0.5 

DISTANCE NORMAL TO BODY/SHOCK-LAYER THICKNESS, y/A 


Figure 43 ELECTRON NUMBER-DENSITY DISTRIBUTION IN SHOCK LAYER 
FOR HEMISPHERE - 20° CONE 
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Figure 44 PLANE-WAVE TRANSMISSION COEFFICIENT AT f = 1.22 
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Figure 47 PLANE-WAVE TRANSMISSION COEFFICIENT AT £ = 50 
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Figure 48 PLANE-WAVE TRANSMISSION 
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COEFFICIENT AT f = 90 






TRAJECTORY POINTS 

FOR WHICH CALCULATIONS 

WERE PERFORMED 

ALTITUDE, ft I VELOCITY, ft/sec 
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TIME FROM 400,000 ft. sec 

Figure 49 TYPICAL SHUTTLE ENTRY TRAJECTORY 












THICKNESS, y/A 

Figure 50 ELECTRON NUMBER-DENSITY 

DISTRIBUTION IN SHOCK-LAYER 
FOR BLUNTED DELTA AT oC = 20° 
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ELECTRON DENSITY, e~/cm 



THICKNESS, y/A 

Figure 51 ELECTRON NUMBER-DENSITY 


DISTRIBUTION IN SHOCK-LAYER 
FOR BLUNTED DELTA AT oc = 20° 
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ELECTRON DENSITY, e cm 



DISTANCE NORMAL TO BODY/SHOCK- LAYER 
THICKNESS, y/A 

Figure 52 ELECTRON NUMBER-DENSITY 

DISTRIBUTION IN SHOCK-LAYER 
FOR BLUNTED DELTA AToc= 20° 
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ELECTRON DENSITY, e“/cm 



THICKNESS, y/A 

Figure 53 ELECTRON NUMBER-DENSITY 


DISTRIBUTION IN SHOCK-LAYER 
FOR BLUNTED DELTA AT od= 20° 



ELECTRON DENSITY, e /cm 



DISTANCE FROM SHOULDER, y, cm 


Fiqure 54 TEMPERATURE IN BASE REGION OF RAM C-IH AS DETERMINED FROM 
SWEPT-VOLTAGE PROBE CHARACTERISTIC 
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ELECTRON TEMPERATURE 
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Figure 55 ELECTRON DENSITY IN BASE REGION OF RAM C-HI AS DETERMINED FROM ELECTRON 
CURRENT REGION OF SWEPT-VOLTAGE PROBE CHARACTERISTIC 




